The 


Canadian Journal 
of 


Chemical Engineering 


formerly 


CANADIAN JOURNAL OF TECHNOLOGY 


CONTENTS 


Theoretical Analogies Between Heat, Mass and Momentum 
Transfer and Modifications for Fluids of High Prandtl A. B. Metzner 
or Schmidt Numbers W. L. Friend 


Flow of Compressible Fluids VW. B. Powley 
High Energy Fuels for Jet Propulsion Engines C. W. Perry 


Mass Transfer in a Bubble Column A. 1. Johnson 
C. W. Bowman 


Recent Developments in the Manufacture of Chlorine 
Dioxide W’. Howard Rapson 


Catalytic Reforming N. J. Enis 


The Manufacture of Chemical Cellulose from Wood VW. Wayman 


The Settling Behavior of Uranium Trioxide-Water Slurries A. W. Boyd 
J. L. Whitton 


Index of Articles — 1958 


Index of Authors — 1958 





Photograph taken at C-I-L Central Research Laboratory, Beloeil, P.Q. 


EXPOSURE:1y second 


This unretouched double-exposure photograph shows 
an individual shock wave from an exploding detonator. 
It was obtained by exposing a film to two successive 
images of the shock wave as revealed by discharging two 
condensers in succession through a spark gap. 

The high intensity light flash for each image was 

of less than one microsecond duration, and the time 
interval between flashes was ten microseconds. 

From such investigations, the efficiency of detonators 
is improved. The end result of research reaches the 
C-I-L customer in the form of improved products, 

and a higher standard of technical service which 
Canada needs for industrial progress. 
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IGHT from the beginning in June 1957, the 

destiny of The Canadian Journal of Chemical 
Engineering has been guided by W. M. Campbell of 
Chalk River, the Journal's editor. 


In October 1958 — nine issues and 53 papers 
later — Editor Campbell wrote “finish” to his duties. 
He had taken on the task of editor on a temporary 
basis to set a standard of quality for the Journal's 
contents; now that the pattern had been well 
established, he asked to be relieved from his edit- 
orial work. 

Having Dr. Campbell as first editor was a most 
fortunate choice for the Institute. Many problems had 
to be solved to keep the issues rolling off the 
presses every two months. Authors had to be per- 
suaded of the merits of publishing in the Journal. 
A list of reviewers had to be set up. And the 
eternal triangle of editor — author — and reviewers 
occasionally created an explosive atmosphere that 
required the firm hand of the editor. How well 
Editor Campbell succeeded in resolving the multi- 
farious problems of the editorial job is shown most 
tangibly by the healthy condition of his “baby” 
today. 

The aim of C.J.Ch.E. is to provide a service to 
all chemical engineers in the country. To do this, 
the editorial board recommended that each issue 
should contain both research and non-research 
papers. The non-research papers could deal with 
subjects such as recent developments in unit oper- 
ations, flow sheets and descriptions of new Cana- 
dian processes, waste treatment and pollution 
problems, equipment experiences, economics and 
operations research. In an editorial in the first 
issue, Dr. Campbell stated that ‘'’C.J.Ch.E. will be 
devoted entirely to chemical engineering and close- 
ly related subjects. For example, applied chemistry 
is a closely related subject if it is slanted towards 
chemical engineering.” 

Aas Editor Campbell predicted in the same issue, 
the problem of securing research papers has proved 
easier than obtaining papers from non-research men. 
However, more papers would be welcomed from 
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the latter group such as plant engineers, design en- 
gineers, operations engineers, etc. 

One of the main problems encountered so far 
with industrial papers is the lack of quantitative 
data in them. Unfortunately, this is a common situa- 
tion in industry where it is not always possible to 
collect all the data desired for a really complete 
paper. 

The 2,700 circulation of C.J.Ch.E. is worldwide. 
About 1,600 members of the Institute receive it 
regularly by choice; and by paid subscription the 
Journal goes to another 130 persons in Canada, 
over 400 in the United States, 75 in the United 
Kingdom, and about 250 in a wide variety of foreign 
countries. 

The October issue included paid advertising for 
the first time. A limited amount of advertising is 
to be carried in the Journal in future issues. These 
advertisements will be particularly suited to an 
engineering audience. Inclusion of advertising is 
expected to help make the Journal's operations self- 
supporting. 

So it is evident that the Journal has come a long 
way under the leadership of Bill Campbell. He 
would, of course, be the first to “play down” his 
role and emphasize that much of the progress of the 
new Journal was made with the help of the authors, 
the reviewers of papers, and the editorial board. 
Without a doubt, all these people made a major 
contribution to the well-being of C.J.Ch.E. However, 
the Journal bears very definitely the authoritative 
stamp of Dr. Campbell's thoroughness, his organ- 
izational ability and his technical and engineering 
talents. 

Fortunately C.J.Ch.E. is not losing his services 
entirely. He has kindly consented to become 
Honorary Editor. In this capacity, his experience 
as pioneer editor of the Journal will prove invalu- 
able to the Institute’s head office staff in Ottawa 
who take over the editorial duties with this issue. 
All future correspondence regarding editorial 
matters should be directed to the Institute’s Ottawa 
address. * 
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Theoretical Analogies Between Heat, 
Mass and Momentum Transfer and 


Modifications for Fluids of High Prandtl 
or Schmidt Numbers’ 


A, B. METZNER? and W. L. FRIEND® 


The presently-available analogies for turbulent 
flow inside smooth round tubes are reviewed in 
detail; an empirical coefficient in the most promising 
equation is evaluated experimentally using the first 
heat transfer data of high accuracy obtained in the 
region of Prandtl numbers above 100. The final 
correlation is based on both heat and mass transfer 
data, and covers Prandtl (or Schmidt) numbers from 
0.46 to 3,000. Up to Prandtl or Schmidt numbers of 
600 the results are correlated with a standard 
deviation of 9.4%. 


— paper is concerned with the rates of heat (or 
mass) transfer to Newtonian fluids flowing turbulently 
inside smooth, round tubes. The familiar empirical equa- 
tions presently used, of the form 


iia SC CN ae Wie iar Seo ce Reba male (1) 


predict various values of the Nusselt number for given 
Reynolds and Prandtl numbers because of disagreement 
as to the numerical values of the constants c,, c, and c,. 
As written, Eq. (1) deals only with the so- called “iso- 
thermal coefficient”, defined as that coefficient which 
would be obtained in the limit of a zero temperature 
difference between the fluid and the tube wall. In practice, 
the required finite temperature difference causes radical 
variations in the physical properties of the fluid. 

A wide variety of methods have also been proposed 
for correlating the differences between such “real” co- 
efficients and the idealized isothermal coefficient (1: 2; 3), 
Again no conclusive choice of a single correction, or 
groups of corrections, is possible on the basis of available 
data. McAdams (2) text reviews this work in some detail 
and concludes that further work is clearly required. 


The development of various theoretical analogies 
between heat or mass transfer rates on the one hand and 
momentum transfer on the other have resulted in a 
similar plurality of design equations which predict appre- 
ciably dissimilar results (+ 5 6 78.9), Recent work (+9) in 
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this area has clearly refuted the concept of a purely laminar 
sublayer next to the wall, but the magnitudes of the 
empirical corrections suggested to account for turbulence 
near the wall depend upon the experimental heat or mass 
transfer rate data upon which they are based. Accordingly 
a clear need exists at the present time for contributions of 
two types: 

(a) A review of the best available analogies between 
heat (or mass) transfer and momentum transfer rates, to 
clearly separate the useful theoretical concepts. 

(b) Precise experimental results are required to fix 
the numerical values of any empirical coefficients which 
may appear in the above relationships. 


The present paper is primarily concerned with the 
theoretical development of an acceptable analogy. A 
companion paper “!) presents new data on heat transfer 
to high Prandtl number fluids and discusses the earlier 
available experimental data on both heat and mass transfer. 


Theory 

Any theoretical approach to the prediction of heat or 
mass transfer rates in turbulent fluid streams must, of 
necessity, proceed from a thorough understanding of the 
local flow phenomena. The laws which govern the turbu- 
lent exchange of momentum are not yet well established, 
however, and an entirely theoretical solution is presently 
impossible. Nev ertheless, much progress has been made 
with simplifying assumptions and hypothetical models of 
turbulent fluctuations, such that it now seems probable 
that the transfer of heat, mass, and momentum may be 
treated with a single theory. 


The equation for momentum transfer may be written 


Ta./p = (Em + ¥) du/dy.... 2c. ce. cccces (2) 
and that for heat transfer: 
q/pCp = — (En + »/Np,) dT/dy............ (3) 


Since heat and mass transfer processes are believed to be 
analogous, Eq. (3) may also be used for mass transfer. In 
this case a concentration gradient dc/dy replaces the 
temperature gradient dT’/dy, the molecular diffusivi ity of 
mass D,, replaces the thermal diffusivi ity v/Np,, a turbulent 
diffusivity E, replaces E,, (although the two are always 
assumed to be equal) and a mass flux replaces the heat 
flux q/pC,. 
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The eddy diffusion coefficients defined by these equa- 
tions, although perhaps of limited physical significance, 
are useful simplifications of the complex process and 
provide a plausible means of relating the various transport 
phenomena. The question which immediately arises is: 
what is the relation between the eddy diffusivities of heat 
and mass (EF, and F,) on the one hand and that of 
momentum transfer (E,,) on the other? The equality of 
these two types of eddy “diffusivities is supported by the 
simpler mechanistic theories of turbulence, such as 
Prandtl’s mixing length hypothesis, ‘'!.!*) and is usually 
assumed by workers in this field. E xperimental confirma- 
tion of this supposition is lacking, however. Measurements 
of temperature and velocity profiles in free turbulent 
wakes (13) indicate a value for E,,/E,, of 2, while the data 
of several investigators for air in pipe flow yield an average 
value of 140. E xperiments with the diffusion of 
gases (15) in the central portions of a rectangular channel 
suggest a value of E,/E,, of approximately 1.6. Yet, sur- 
prisingly enough, the best agreement between theory and 
heat transfer data is obtained when E,,/E,, is taken as 
unity (7), 

Reichardt offers a possible explanation of these seem- 
ingly incompatible observations by noting that the experi- 
mental evidence that F,,/E,, is smaller in friction layers 
than in free turbulence permits the tentative conclusion 
that the reduction of this ratio is a result of the influence 
of the wall. It follows that E,,/E,, decreases with decreas- 
ing distance from the wall. Reichardt suggests that E,,/E,, 
is unity in the immediate proximity of the wall and tends 
to approach the free turbulence value of 2 at larger 
distances from the wall. 

Reliable determination of the diffusivity ratio distri- 
bution across a tube requires extremely precise measure- 
ment of the temperature and velocity fields, and only a 
few notable investigations have been performed. Isakoff 
and Drew ‘!%), in a study of heat transfer to mercury 
found that E, /E ‘m did increase with distance from the 
wall in the immediate vicinity thereof, but went through 
a maximum at a value of y/R of approximately 0.20. The 
available data for air are summarized in a recent paper by 
Sleicher (4), 

Before the basic transport equations may be employed 
to compute rates of heat transfer from flow measurements 
several assumptions must be introduced: fluid properties 
must be assumed constant, or specified as some arbitrary 
function of temperature; the local heat-flux must be 
described as a function of distance from the wall; a 
relation between the eddy diffusivities for heat and 
momentum must be proposed; and, finally, the eddy 
diffusivity for momentum must be expressed analy tically 
in terms of the relevant variables. 

Restricting the theoretical developments to the iso- 
thermal heat transfer coefficient is not an especially 
important limitation at this time, particularly in view of 
the uncertainty introduced in other assumptions. In fact, 
it is more likely to be preferred for the sake of simplicity. 
Much the same is true of variations in the heat transfer 
coefficient with tube length, L/D. Since this factor has, 
however, been clearly defined ‘'”) to be absent from tubes 
having an L/D ratio greater than about 20, elimination of 
this complication is of even less importance. The local 
heat-flux may be taken with a high degree of accuracy 
as linear in radius: : 
G/Go°@ 1 = Y/R wv ccccvaecks ....(4) 


This approximation depends upon the assumption of uni- 
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form temperature and velocity profiles, and the error 
introduced, as estimated by Reichardt (7) is small except 
for extremely low values of the Prandtl number. A similar 
equation for the shear distribution may be derived 
rigorously (#18), 


Blip ERP Ree sie dinia hid sista nic ois 


By far the most critical of the assumptions is intro- 
duced in the analytic description of the eddy diffusivity 
variation. In the absence of an acceptable theory, Eq. (2) 
and (5) must be employed simultaneously to compute E,, 
from experimental measurements of the velocity distribu- 
tion. This procedure, which yields reliable results in the 
turbulent core of the flow, fails badly in the important 
region near the wall where the necessary precise measure- 
ment of point velocities is extremely difficult, if not 
impossible. Consequently, in this region the nature and 
magnitude of the turbulence must be specified entirely 
arbitrarily. Since the eddy diffusivity must be zero 
exactly at the wall it must be very small in the immediate 
vicinity thereof. However, it may not be negligible, 
especially for moderate and high values of the Prandtl 
number. This may be illustrated by considering Eq. (2) 
and (3) for a fluid with a Prandtl number of 100. Even 
if the turbulent friction is only one hundredth of the 
molecular (E,,/v = 0.01) the turbulent and molecular 
heat-fluxes will be of the same order of magnitude: 


Eh 0.01 
vy/Np, _-v/100 


The assumption of equality of the eddy diffusivities of 
heat and momentum in this calculation, is seen not to be 
critical in arriving at the conclusion that the turbulent 
heat flux is very important. 

Several authors have proposed solutions based on a 
particular assumption for the eddy diffusivity variation, 
but perhaps the most significant contribution to this field 
has been made by Reichardt ‘7) who succeeded in deduc- 
ing the general form of the solution prior to the inclusion 
of any hypothetical suppositions regarding the nature of 
the flow. This dev elopment, free of any major assump- 
tions, progresses the art one step further than the con- 
ventional dimensional approach, revealing not only the 
proper grouping of variables, but the formal relations 
involved. A resume of his derivation is given in reference 
(1%), In view of its importance to the present problem and 
the comparative unavailability of references “) and (, 
the basic steps are outlined in the following paragraphs. 


Upon introducing a dimensionless temperature differ- 
5 ae 5 : . 
ence, 6 = -.; =. in the equations for the point and 
r¢ — To 
total heat-fluxes one obtains the following equation for 


the dimensionless temperature gradient: 


dé q 1 dé 
; = ; : pee WE oS leeds (6) 
d (y/R) do \1 + N*p En/v d (y/R)/, 
ee En 
= Nw (5°) 


signifies a turbulent Prandtl number. 


Here N*>, 


Similarly after introducing the dimensionless velocity 
¢@ = u/U in the equations for the fluid shear, there follows 
for the velocity gradient: 


dp of I dd 


d(y/R) ~ AT + Ewe) aq BJ 
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Dividing Eq. (6) by (7) and nothing that both q and + 
may be taken as linear in radius (Eq. (4) and (5)) one 
obtains: 

dé 1 + En/v dé 
dp 1+N*p, En/? dg), 


Eq. (8) is the basic differential equation governing the 
relationship between the temperature and velocity fields. 
For the computation of the heat transfer, the Stanton 
number is related to (d¢/d@), by: 


N a dé f/2 Pm/ Om 
\dp), Nex 


The quantities $, and 6, represent the ratios of the mean 
to the maximum velocities and temperature differences, 
respectiv ely. 


The dimensionless derivative (d@/d¢), may be ob- 
tained by integration of Eq. (8) from the wall to the 
centre-line of the tube: 


1 1 
-, = (99) J f(_1+ Eni l 
foo 1 = (24) { [222.0 
°o o 
1+ En/v 
er i/ | (; + Ne wr). 


oO 


By algebraic manipulation this result is simplified to: 


(3) a ca 
1 a 1 
dd : ‘ 


3 ¢ 
Np,* — 1 eae 
1 + (Nz | + Ne.* E,/p 


oO 


Anticipating the use of the universal velocity distribu- 
tion equation, it is desirable to change the variables of 
integration from ¢ to ut = u/u’: 


dé Np,* 
aaee YN *digk’ caresses —....(13 
(i), ue = 


dut* 
Np,* — 1) u*/U 
7 ; Te | aes En/v 


o 


The ratio u*/U is given in more familiar terms by: 


u*/U = om Jte/v Dip IDs ences 1s 
p 


Combining this result with Eq. (9) for the Stanton 


number: 


f/2 7 (Ner*/Ner) 


U/u* 


dut 
1 + Np,* 


o 


1 + (Np;* = 1) dm 


U/u* 
ie Si dut 5 ‘eit ie alii 
The integral, 1+ Ne E./p in which the eddy 
Oo 
diffusivity variation has been isolated is denoted by 
Reichardt as “b”. 
Thus, 


= (f/2) (dm, Om) (E n/ em) 
Nst = ae 


1 + (Np;* ny 1) dm 2 (b) 
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One may note that the integrand of the defining integral 
for b is identical to the point value of the fractional heat 
flux which is due to pure molecular conduction: 


Qmo 1 1 1 a 
= — = — = - ——...(17) 
do 1 + (q:/qmo) Cp9En 1 + Ner* Em/v 
1+p k 


Finally, introducing the assumption that E,, = E, the 
desired results are summarized by means of the equations: 


; : (f/ 2) (dm /Om) : 
1 + (Npr — 1) (om) Jin (b) 


U/u* U/u* 


dut* ( Qmo 
b = - = ——— = 
1 + N Pr Ba v | q 


Oo Qo 


Nst _ 


Engineering simplification of the theory 

The final results obtained by the methods of the 
previous section, Eq. (18) and (19), are still quite com- 
plex, particularly from the point of view of ready engin- 
eering application. Fortunately, the solution can be 
simplified further for most applications with approximate 
representations and by excluding the region of very low 
Prandtl number Np,<<1, (ie., heat transfer to liquid 
metals) from consideration. To this end, in this section 
the variation of the individual parameters will be examined 
in detail. 

The quantities ¢,, and 6, are terms of secondary 
importance, both being close to unity. The ratio of the 
mean to maximum velocities, },, increases slightly with 
the Reynolds number, varying from 0.78 to 0.84 over the 
usual range encountered, 10*<N,p,<10®, ‘?). The ana- 
logous temperature difference parameter, 6,,, increases 
with both increasing Reynolds number and Prandtl 
number, approaching a value of unity asy mptotically for 
high Prandtl numbers. For unit Prandtl number it is 
closely equal to ¢,, and, on the basis of experimental 
temperature profiles for water (Np, = 4.0) reported on 
p- 208 of reference ‘*), the rate of approach to unity with 
increasing Prandtl number appears to be rapid. For 
Prandtl numbers much less than unity 6,, changes rapidly 
with both the Reynolds and Prandtl groups and can no 
longer be treated as a second order parameter. 


Thus, with the exception of the region of low Prandtl 
numbers, the major problem of the turbulent heat transfer 
has been resolved to one of evaluating the Reichardt b 
function of Eq. (19). The integrand of the defining equa- 
tion for b represents a ratio which will have a maximum 
value of unity exactly at the wall and will decrease with 
increasing distance from the wall, approaching zero in 
the turbulent core at a rate determined by the Prandtl 
number. Again, with the exclusion of the low Prandtl 
number region, the ratio will be of significant magnitude 
only in the immediate vicinity of the wall, and the b 
function will be determined exclusiv ely by conditions near 
the wall. In this region the turbulent exchange, expressed 
as the eddy diffusivity, — only upon the dimension- 
less distance from the w all, yt, and is independent of the 
bulkstream conditions or Rey nolds number (7), Alter- 
nately, one may arrive at the same conclusion by notin 
that the eddy diffusivity may be expressed solely as a 
function of ut, the variable of integration in Eq. (13), 
as ut+ and y+ are uniquely related over the entire cross- 
section of the tube. These observations permit the impor- 
tant conclusion that b is a single-valued function of the 
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- 


Prandtl number and is independent of the bulk stream 
flow or Reynolds number. 


The necessary conditions for this conclusion con- 
cerning the variation of b will generally not be met for 
low values of the Prandtl and Reynolds numbers, and it 
is important to formulate numerically approximate limits 
of applicability. For this purpose assume that the universal 
velocity distribution is given by 


De ales tk eames Boe ale » 20) 


where g is an arbitrary differentiable function of the 


argument y+. Then, 
dut r( 
Ay"? 
dy* 
and 
du Toke du* Toke 
= — = PREY, ow cio tdannaeeee 
dy py dy* pv 


Combining Equations (2), (5) and (21) gives: 
En (1—y/R) 


v eAy" 


Substitution of this result into Eq. (19) gives upon 
defining* the “wall region” by 0<y/R<0.10, and assuming 
that 1 — y/R is equal to unity over this region: 

y/R = 0.10 YRt 

g’ (y*) dy? g’ (y+) dy 

b= + —— (23) 


O 1 y/R = 0.10 _ 1—v/R 
1 + Np, — 1 1 + Np, “$e — 1 
g’ (yt) g’ (yt) 


where y,* is the value of y+ at the centerline, y = R. 


The second integral in Eq. (23) introduces the 
Reynolds’ number via the tube radius, R, and the upper 
limit y,+. Therefore, the function b itself will be inde- 
pendent of the Reynolds number only if the second 
integral is negligible or zero. This condition may be 
formulated approximately by restricting the intergrand 
to negligible magnitude at the boundary of the wall region, 
y/R = 0.10. 


g’ (y*) 


; 1  O, say 10-2... .(24) 
a lee = 4 
g ty") y/R = 0.10 


The maximum values of the integrand is unity exactly at 
the wall, and therefore, the limit of 10-2 selected is con- 


servative. 


Differentiating the logarithmic velocity law ‘) to esti- 
mate g’(y+) gives: 
dut 
= gy’ (y+) = 2.5/y*.. ; .(25) 
dy* . : 
Substituting this result into Eq. (24) and neglecting 1 in 
comparison to y + /2.5 and 1 in comparison to Np,(y +/2.5) 


one obtains: 
. fy 
wT :) | > 100....... (26) 
2.5 y/R = 0.10 


*At the limiting, lowest value of Nre of interest, 10*, this definition vields 
in terms of y+, 0 < y+ < 32, which corresponds to the regions usually referred 
to as the laminar and buffer layers. At Nre = 105 it implies O < yt < 240, 
which is more conservative than necessary. 
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From the definition of y+: 


Y  |Tofe 
(yt)y/R = 0,10 = 
y y V ? 


_ 0.10R f V2 
= : : 
= 0.05 Nre y' Be 2 cotta ale (27) 


The final condition is then: 
(Npr) (Nre)? (f£) See Be a os ey cols (28) 


For high Prandtl number, this condition will be met 
for all turbulent Reynolds numbers, but for lower values 
of Np, the condition requires slightly higher Reynolds 
numbers. This is not an especially limited restriction, 
however, as operation with high Reynolds numbers is not 
only feasible, but generally economical for non- viscous, 
low Prandtl number fluids. For a lower Reynolds number 
of 10?: 

ig: cc dc ee (29) 


Reichardt has performed a similar analysis in which 
he imposed the restriction that qmo/q fall to the value 
of 0.50 for y/R<0.05. His result with this less severe 
restriction was: 

Ne Nae > 2500......... =i) eo er 


Or, for a lower Reynolds number of 10+: 
Np, > 0.25 


Semi-theoretical solutions 


Since Reynolds first postulated (2°) that heat and 
momentum are transferred by identical mechanisms, much 
attention has been focused on the analogy, and several 
authors have proposed more refined solutions to the prob- 
lem; these have been reviewed in detail in reference ") in 
terms of the general considerations of the prior sections. 
To conserve space, only the three best types of approach 
will be mentioned here and the reader is referred to the 
literature ‘4 2,12,21) for more complete reviews. 


Von Karman ‘®) subdivided the flow into three zones: 
the laminar, buffer, and outer turbuelnt layers. Utilizing 
the velocity distribution data of Nikuradse, expressions 
were developed between u+ and y+ in each of the flow 
zones, and the eddy diffusivity for momentum was deter- 
mined by differentiation. The ratio of eddy diffusivities 
in the heat transfer equation, ¢,,, and @,, were each taken 
as unity. von Karman’s result for b is: 

_| ‘ In (1 + 5/6,[Nvr — 1])| 


b = 551 


Dak gninte a 31) 
| Np, — 1 ( 


Matinelli ‘®) extended the approach of von Karman by 
including the effect of the main-stream conditions, or 
Reynolds number, on b. His result, as expected, does not 
differ appreciably from von Karman’s except for values 
of the Prandtl number much less than unity. 


The equation of von Karman fits the available data 
quite well up to Prandtl numbers of about ten; for higher 
Prandtl number it underestimates the heat transfer. This 
is a consequence inherent in the model of an entirely 
laminar layer adjacent to the wall. For high Prand¢l 
number all of the temperature gradient is found within 
this layer, in which by assumption, only molecular transfer 
can occur. For this limiting case von Karman’s equation 
reduces to: 


Nse 
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The identical result is obtained by the crude film theory 
when the film thickness is obtained from the usual 


- 


definition of the so-called laminary sub-layer, y+ = 5. 


The next logical extension of the theory is, of course, 
to rectify the model of an entirely laminar zone by 
developing a functional relation for the eddy diffusivity 
which vanishes exactly at the wall rather than at the 
boundary of a superficial laminar layer. As previously dis- 
cussed, the presence of a small amount of turbulence in 
this region would not significantly affect the velocity 
rofile, but would markedly improve the heat transfer at 
high values of the Prandtl number. Unfortunately, how- 
ever, at the present time the nature and magnitude of this 
critical turbulence cannot be deduced from basic flow 
measurements as in the turbulent core. In fact, the only 
quantitative experimental evidence for its existance is the 
apparent disagreement of the available data for heat and 
mass transfer rates to fluids of high Prandtl or Schmidt 
numbers with computations based on a purely laminar 
layer. Several authors have proposed solutions based upon 
an arbitrary designation of the eddy diffusivity variation 
in the nearly laminar region near the wall, but these 
authors’ results can hardly be regarded as anything more 
than semi-theoretical correlations of data achieved by 
arbitrary assumptions. The analogies of Deissler ‘*); Lin, 
Moulton and Putnam (5); Rannie (*); and Reichardt “) are 
all typical of this technique. 


Since meaningful flow measurements extremely close 
to the wall are not yet available, the only approach 
possible—other than use of arbitrary assumptions as men- 
tioned above—is to use the heat transfer data to formulate 
the Reichardt b function experimentally. By excluding the 
region of the low Prandtl number from consideration, b 
is determined exclusively the Prandtl number. Therefore, 
analysis of the data at one Reynolds number suffices. 
Simplifying assumptions applied to enable use of the Eq. 
(18) are as follows: 

(a) The diffusivity ratio 


taken as unity. 


E,,/Em, and 6,, were both 


(b) The mean to maximum velocity ratio was assumed 
constant, at an average Vv value of 1/1.2, which is approxi- 
mately representative of a large range of Reynolds num- 
bers. Inserting these assumptions into Eq. (18) gives 

f/2 
1.20 + (Nw — 1)b “ao 


where b is to be determined empirically. Eq. (33) will be 
valid for fluids with Prandtl numbers greater than 0.25- 
0.60, as discussed earlier. 


Ns = (33) 


Experimental results: evaluation of the b function 
and correlation of data 


The correlation of experimental data is described 
detail in the companion paper “'*), However, for com- 
pleteness the salient results are presented here as well. 
The final correlation was based on the following data: 

1. All experimental heat transfer results of high 
accuracy, excluding only liquid metals, which were avail- 
able in the literature were utilized. These data covered a 
200-fold range of Prandtl numbers from 0.46 to 95. 


2. To extend the correlation into the region of high 
Prandtl numbers of primary interest here, experimental 
data on viscous syrups and molasses solutions were taken 
over about a 10-fold Prandtl number range from 55 to 590. 

3. Approximate mass transfer data, covering a range 
of Schmidt numbers from 1,000 to 3,000, and also available 
in the literature, enabled the extension of the correlation 
—at least approximately—to Prandtl or Schmidt numbers 
of 3,000. 
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The value obtained for b was 
ee eT (34) 


Insertion of Eq. (34) into (33) gives the final corre- 
lation: 
{/2 
Ns = — : eat re (35) 
1.20 + 11.8 Wf/2 (Ne — 1) (Ne)-¥3 


The maximum deviation of any of the heat transfer 
results from the correlating equation was 22%; the 
standard deviation was 9.4%. Inclusion of the Sieder- 
Tate (3) viscosity correction to reduce errors due to use 
of finite temperature differences in the experimental 
results taken from the literature would have reduced 
these deviations to 15 and 6.5%, respectively. Thus, a 
highly accurate design procedure has been developed 
over a 1300-fold range of Prandtl or Schmidt numbers to 
600; it has been extended at least approximately to en- 
compass the 6500-fold range from 0.46 to 3,000. 

A design chart based on the final correlation (Eq. 35) 
is given in Figure For mass transfer the coordinates 














Samana aa 


Figure 1 





become Ngy (k,/V) and Nge (u/pD,) ) instead of Ng, and 
Np,, respectively. The line shown at Nag, = 10* is the 
actual line through the experimental data; the curves at 
other Reynolds numbers, (at Prandtl numbers above 100) 
represent extrapolations out of the region in which the 
experimental data were obtained. This introduces the 
possibility of uncertainties in use of the correlation at the 
highest Reynolds numbers when the Prandtl or Schmidt 
number is also high. While this limits the prediction 
of mass transfer rates in the case of heat transfer 
Reynolds numbers of 3 x 10‘ and higher are extremely 
difficult to obtain if the Prandtl number is above a few 
hundred because of the highly viscous fluid systems 
involved: pressure gradients of the order of 10 p.s.i. or 
more per foot of tube length would be required with 
fluids in common usage. 


Conclusions and design recommendations 


An analogy, based upon earlier work by Reichardt “), 
between heat, mass and momentum transfer is presented 
which allows for influences due to slight turbulence in 
the regions close to the wall. The analogy is formulated 
in general terms, without the introduction of any arbitrary 
assumptions as to the nature of this turbulence. By 
excluding only heat transfer to liquid metals, the theo- 
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retical equation could be greatly simplified. All other 
available heat and mass transfer data were used to evaluate 
the single empirical coefficient in the equation. The final 
recommended design equation for heat transfer is: 


f/2 
Nst — 
1.20 + 11.8 Aff/2 (Nee — 1) (Ne)-¥8 


For mass transfer, Ng, replaces Ns, and Ng, replaces Np,. 


The above equation was based on heat transfer data 
on fluids having Prandtl numbers between 0.46 and 590. 
Approximate extension to Prandtl or Schmidt numbers of 
3,000 was achieved through use of mass transfer data. A 
design chart based on this equation is given in Figure 1. 


The above predictions of the “isothermal” heat trans- 
fer coefficient must be suitable modified for use under 
conditions of high temperature difference between the 
tube wall and the fluid. If the fluid viscosity is the only 
variable which changes rapidly with temperature, use of 
the Sieder-Tate ‘) (m,/u,)°!* or the Kreith-Summer- 
field “') (40,/uy)°!° corrections is recommended. 
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Nomenclature 


Any consistent set of units may be employed. Those listed 
are merely illustrative. 
U/u* 
= Reichardt b function, b = | — 
J q 
Oo 
mentally by 11.8 (Np,)~'/3, dimensionless 
C1,C2,c; = Empirical constants, dimensionless 
Specific heat, B.t.u./lb. mass °F. 
Differential operator 
Tube diameter, ft. 
Liquid molecular diffusion coefficient, ft.2/hr. 
Eddy diffusivity for mass diffusion, ft.2/hr. 
Eddy diffusivity for heat, ft.?/hr. 
Eddy diffusivity for momentum, ft.2/hr. 
Fanning friction factor, dimensionless 
Functional representation of universal velocity distri- 
bution, ut = g (yt) 
Conversion factor 32.17 (Ib.m) (ft.)/(Ib.#) (sec.?) 
Mass velocity, lb.m/(hr.) (ft.?) 
Film coefficient of heat transfer, B.t.u./(hr.) (ft.2) (°F.) 
Thermal conductivity, B.t.u./(hr.) (ft.) (°F.) 
Liquid phase mass transfer coefficient, moles/(sec.) 
(ft.2) (moles)/ft.3 
Tube length, ft. 
Heat flux, B.t.u./(hr.) (ft.?) 
Tube radius, ft. 
Temperature, °F. 
Time average point velocity in x direction, ft./sec. 


) du*, given experi- 


= Friction velocity, Yg.7/p, ft./sec. 


ut = Generalized velocity parameter, u/u*, dimensionless 
U Maximum linear velocity at centerline, ft./sec. 

V Mean linear flow velocity, ft./sec. 

y+ Generalized distance parameter, yu*/v, dimensionless 
yrt = Value of yt at axis of tube 

x, y, Z = Cartesian distance coordinates, ft. 


Greek letters 


Dimensionless temperature difference T — T./Tg¢ — T, 
= Ratio of mean to maximum temperature differences 
Viscosity, lb.m/ft. sec. or lb.x/(ft.) (hr.) 
Kinematic viscosity, ft.?/hr. 
= Fluid density, lb.s/ft.* 
Shear stress, lb.¢ /ft.? 
Dimensionless velocity, u/U 
dm Ratio of mean to maximum velocities, V/U 


Subscripts 


b = Bulk average 
¢ Certerline 
mo Molecular 

oO Wall surface 
w Wall 


Dimensionless Numbers 


Nnu = (hD/k) Nusselt number 

Npr = (C,u/k) Prandtl number 

Np-* = (Np-) (En/Em) Reichardt turbulent Prandtl number 
Nre = (DG/pu) Reynolds number 

Nsc = (u/oD_.) Schmidt number 

Nsh = (k,/V) Sherwood number 

Nst (h/C,G) Stanton number 
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Flow of Compressible Fluids’ 


M. B. POWLEY- 


Previously derived equations describing compres- 
sible flow have been rearranged in terms of known 
downstream conditions to facilitate determination of 
unknown upstream conditions. A 650 I.B.M. com- 
putor was used to solve the simultaneous equations 

/ numerically to permit construction of convenient 
design graphs in parametric form. Use of these graphs 
will eliminate a tedious trial and error computation 
formerly necessary with presently available methods 
when the upstream pressure and temperature are to 
be determined in gas flow problems. Application to 
design is illustrated by a worked sample. 


HEN designing pipe lines for handling compressible 

fluids, engineers often follow this “rule of thumb”: 
If the pressure drop is less than 10% of the upstream 
pressure use an average density for the fluid and then 
continue calculations with incompressible flow equations. 
If the pressure drop is greater than 10% of the upstream 
pressure, then compressible flow design methods should 
be used.” It is this latter case that gives rise to consterna- 
tion among engineers, especially those that have to do 
such work occasionally. It is generally believed that the 
equations describing compressible flow are complicated 
and unwieldy, and are best avoided if possible. One of 
the aims of this report is to dispel this idea and to show 
how exact compressible flow calculations can be per- 
formed in a simple, easy manner. 

The condition of high pressure drop in compressible 
flow frequently occurs in venting systems, vacuum distil- 
lation equipment and long pipelines. In many cas:s, the 
pressure drop is critical, thereby requiring accurate 
analysis and design. 

Theoretically, the flow can be either isothermal or 
adiabatic, depending upon heat transfer through the pipe 
wall. In the usual case that occurs in chemical plants and 
refineries, the piping is insulated, so that the heat transfer 
is nil and the flow is essentially adiabatic. Unfortunately, 
the equations describing adiabatic flow are complex, w hile 
those for isothermal flow, which seldom occurs, are 
relatively simple. It is this fact that has largely deterred 
process engineers from using correct compressible flow 
calculation methods. 


1Manuscript received September 3, 1958. 
2Textile Fibres Department, Du Pont Co. 
Maitland, Ont. 

Contribution from Du Pont Co. of Canada (1956) Limited, Maitland, Ont. 
Based on a paper presented at the Joint A.I.Ch.E.-C.1.C. Chemical 
Engineering Conference, Montreal, Que., April 20-23, 1958. 


of Canada (1956) Limited, 
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In 1943, C. E. Lapple “) derived compressible flow 
equations, which he presented graphically in parametric 
form. To use Lapple’s design graphs, it is necessary to 
know the upstream conditions of pressure, temperature, 
Viscosity, and ratio of specific heats at constant pressure 
and volume. Then the pressure drop in a pipe line can be 
determined for a given flow rate, or vice versa. If the 
downstream conditions are known and it is required to 
determine the upstream conditions, a trial and error pro- 
cedure is necessary with these graphs. As this latter case 
is quite common, it was believed that a set of design charts 
based on known downstream conditions would prove use- 
ful to designers. This work has been completed and is 
reported in this paper. It should be noted that the design 
methods developed here do not supplant those of Lapple, 
but rather complement them. 


Reason for the work 


This work was initiated during a study of the venting 
capacity of one of the process buildings at the Du Pont 
Co. of Canada Nylon Intermediates Plant. A trial and 
error procedure with Lapple’s graphs was first adopted, 
and then abandoned after it was found to be hopelessly 
inadequate for the complex, highly branched, venting 
system of the plant. Lapple’s equation for isothermal com- 
pressible flow was rearranged in terms of known down- 
stream conditions to permit construction of new design 
graphs. This allowed the use of atmospheric pressure as a 
starting point for calculation of upstream pressures at the 
various pipe junctions and in the many process vessels. 


Later, the isothermal design graphs were found to be 
= in design and analysis “of vacuum distillation units. 
Vapor pipes and condensers, operating at pressures. down 
to 5 mm. Hg. absolute, were both designed and checked 
by means of the isothermal flow graphs. 


At very low pressures, it is essential that vapor pipes 
and condensers be designed for very low pressure drops. 
This requires accurate calculations, therefore, errors in- 
herent in the method have to be reduced to a minimum. 
As the systems studied in this work operated under 
adiabatic conditions, assumption of isothermal flow did 
involve some error in the calculations. For this reason, it 
was decided to complete the above study by including 
adiabatic flow of compressible fluids. 


Previous methods of calculation 





Lapple ‘!) made a very complete summary of available 
methods of calculation. He discussed experimental data 
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that proved the friction factor is independent of Mach 
number up to acoustic velocities, and that accepted 
friction factor-Reynolds number graphs can be used for 
both compressible and incompressible fluids under all 
conditions of flow. 


Design equations and graphical methods advocated in 
older textbooks generally were based upon simplifying 
assumptions that incurred large errors at high velocities. 
A reliable method was developed by Stodola °), which 
required plotting of “Fanno lines” on a temperature- 
entropy diagram. This has not gained wide acceptance 
because of the tedious nature of the calculations. 

The major advance in the evolution of compressible 
flow calculation methods was made by Lapple. His design 
charts have been published by Perry “), Brown (#4), and 
others, and are now considered standard chemical engin- 
eering design methods. Use of Lapple’s graphs will not 
be explained in detail here, as this has been well covered 
in the reference works mentioned above. 


Proposed method 

The basic flow equations developed by Lapple have 
been rearranged in terms of known downstream fluid 
conditions of pressure, temperature, and specific volume. 
This allows calculation of upstream conditions, and hence 
line pressure drop, if the mass flow rate is known. In the 
case where the pressure drop is known, and the flow rate 
is to be determined, either method may be used. 

There is one major difference between the two 
methods. In Lapple’s method, the entrance effect of a 
frictionless nozzle is included in the design graphs. In 
the proposed method, the calculations concern the flow 
between two points in a length of pipe. Entrance effects 
have to be allowed for separately. 

In the design graphs presented here, P,/P, is plotted 
against G/G, for various values of N. There are three 
graphs, Figures 1, 2, and 3, for k = 1.0, 1.4, and 1.8 
respectively. All of the graphs pertain to adiabatic flow 
conditions, but the graph based on k = 1.0 may be used, 
without error, for isothermal flows. 

The quantity, G,, which is the sonic mass velocity, 
is us.d to characterize the fluid at the downstream 
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Figure 1—Compressible flow graph for k = 1.0. 


242 


eh 


G=cp 


= Pp vi 
Pv 
SS lee 
a >} kgMw 
{ RT 
S kgP 
7 


P (mm Hg) 
3.34 





1 \WAWNA 





M 


AVI 


// 
X 


D 
V/ 


D 


KA 


‘| 


4 


W) 
V 


J 


Oh 


Jy 


TINNS 
NS 


Y 


Ws 
7 


AX 


VY, 


ce 
dp 


7 
a? 


Z| 


we 


TAWS 


CAN 


> 


Ch 


4 

















100 N 
100 
a2 eI 
ou ‘ai 
Ta 
— lass joss 
qT 
0 ol 02 aS OF 


FA 


\NNS 


“ hee IN 











ol 
Tz 
T | 099 098 a9 


° ol a2 0.3 oO. 








i 
7 


wn 


MENA 


ny 


“ 


HLT VV V7 
AE 





reference point in the pipe. The relationship between G, 
and physical properties of the 
following expressions: 


fluid are given by the 


Ib./(sec.) (sq. ft.) 


U1 
XL 


Vy 
mM 
Gi 
VX 


LW/) 
YA 


M7) 


rhs 
th 


<mIW/) 
il tH 


°} 
NI 
2 
° 
se 
oO 


A 
-—, 
ie 
p< ol 
oz. ON 
04 INF 


no 
DKA 


iL. 
LAM 


Y 


Ly] 
Ly 
ttt 


i 





ML, 
KA) 


tnt 


fH 
Wi 





4 as 0.6 
MACH. No. G/Gc 


2) 
° 
e 
° 
5 


Figure 3—Compressible flow graph for k = 1.8. 
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Other variables, P,/P, and N represent the ratio of 
downstream to upstream pressure and the number of 
velocity heads pressure drop, respectively. The latter term 
js evaluated by summing velocity head loss contributions 


fL 


from friction , entrance effects, and pipe fittings, etc. 


Included on the graphs are lines of constant temperature 
ratio, so that the upstream temperature may be found. 

If the mass flow rate in a given pipe is known, the 
Reynolds number, friction factor, N, and G,, can be 
calculated. Use of these values with the appropriate graph 
will allow determination of the upstream pressure and 
temperature. If the value of k for the flowing fluid does 
not coincide with a value on the design graphs, it is 
permissible to interpolate between graphs. 


The design graphs 

The compressible flow equations upon which the 
design graphs are based are derived upon the following 
assumptions: 

1. Bernoulli’s Theorem 

2. The Ideal Gas Law 

3. A constant friction factor 

4. Horizontal piping 

5. Adiabatic flow 


Two equations are required to describe the flow rela- 
tionships. These comprise an equation of state and a 
mechanical energy balance equation. The adiabatic equa- 
tion, PvK = constant, applies only to a non-flow frictionless 
expansion or compression and does not take into con- 
sideration the changing kinetic energy of a flowing com- 
pressible fluid. An exact equation of state can be found 
by combining a thermodynamic relationship with an 
overall energy balance equation. The derivation of the 
two equations, upon which the design graphs are based, 
is given in the appendix. The equation of state is: 


mente[te-Gy) 


P2 “v1 


The flow equation is: 


2kN 2+ (k — 1) M2? v:\' cai vs . 
> = =— 31 
k+1 (k + 1) M2? Vi "AV, 


Elimination of the specific volume ratio, v,/v,, from 
these two equations would provide the desired relation- 
ship between pressures, Mach number and the friction 
term. While theoretically possible, this is not practicable 
because of the complexity of the equations. However, a 


N and k 


: ; ; Vy 
numerical solution for M, in terms of —*, P,/P,, 


1 
has been carried out with the aid of an IBM 650 computer. 
In addition, the ratio T,/T, has been determined from 


Pave T. 


the Ideal Gas Law, =—, 
Pivi lr; 

General discussion 

In Lapple’s paper, a detailed discussion is given on the 
subjects: 

(a) Comparison of isothermal and adiabatic flow 

(b) Approximate equations for maximum discharge 

(c) Flow of non-ideal gases 

(d) Effect of fittings and valves 

(e) Effect of inlet shape 

(f) Velocity distribution in the pipe 

(g) Comparison of calculated and experimental flow 

data. 
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A very brief summary of Lapple’s comments will be 
outlined here; reference to the original article should be 
made by interested persons. 

The maximum difference in discharge rates between 
isothermal and adiabatic flow will never exceed 20%. 
This maximum difference occurs for short frictionless 
pipes. With the long pipes (N = 50) there is no appre- 
ciable difference between the two types of flow. 

Approximate equations have been developed to des- 
cribe maximum discharge rates for long pipes (N = 50) 
based on known upstream conditions. These are: 


= g Po 
, N Vo 
“it g Mw 
R ToN 


Lapple estimates maximum errors of 6% and 2% for 
values of N equal to 50 and 200 respectively, with these 
equations. 

If non-ideal gases are encountered, the design graphs 
may be used if known acoustic velocities are used in 
calculating G,. An allowance for this case can also be 
made by multiplyi ing the gas constant R by the compres- 
sibility factor. 

The effect of fittings and values in a pipeline can be 
taken into account by adding the number of velocity 


ei fL. 
heads pressure loss to the friction factor term, D An 


alternate process is to add an equivalent length of straight 
pipe to the actual length before calculating N = The 
necessary data can be found in standard chemical engin- 
eering text books. 

Allowance for entrance and exit effects can be made 
in the same way as for incompressible flow. 

The flow equations described here are based on a 
uniform velocity distribution. As this is not quite true, 
a small error should be introduced in the kinetic energy 
part of the derivation of the equations. Experimental data 
quoted by Lapple shows good agreement with values 
calculated by his method. Therefore, the effect of non- 
uniform velocity can be ignored and the design graphs 
can be used with confidence. 


Application to design 


Use of the design graphs should be self-evident. Lapple 
gives worked numerical examples to which reference 
can be made. The only difference between the two sets of 
graphs is the point of reference, i.e., known downstream 
vs. known upstream conditions. To demonstrate the use of 
the graphs, a typical plant problem is worked out in 
the accompanying example. 

Example problem 

The problem is to determine the pressure in a vacuum 
evaporator for a flow of 4,000 lb./hr. of an organic vapor 
through a 100 ft. vapor pipe to a condenser head which 
is maintained at 5 mm. Hg. absolute. An entrainment 
separator is included in the pipe run. The pipe is 20 in. 
sch. 40 and includes two long-sweep elbows. Because the 
ratio of specific heats is nearly unity, the flow can be 
taken as isothermal. 


Additional data 

Molecular weight —100 

Viscosity —0.014 centipoises 
Ratio of specific heats—1.03 

Absolute temperature —450°K. 

Pipe I.D. —1.568 ft. 

Pipe flow area —1.931 sq. ft. 
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Calculations 


M A G 4000 lb. 

. Mass flow rate, ex 

3600 (1.931) (sec. (sq. ft.) 
= 0.576 


2. Acoustic mass flow rate P mm Hg. 


[| k Mw 
at 5 mm Hg. G, 3.34 


VT (CK) 
5.0 /1.03 (100) 
3.34 Yo 450 


eo Ib. 
= 0.716 
(sec.) (sq. ft.) 


. aie 0.576 
3. Ratio, G/G, ee 
= ().804 
4. Friction factor, f. 
DG 
Reynolds No. ie tas 
. be 
1.568 (0.576) (1488) 
2 0.014 
= 96x 10! 
Therefore, friction factor, 
f. = 0.020 


{L 0.02 (100) 
D 1.568 


5. Velocity heads friction term, N. 
Entrance loss = 1.5 vel. heads 
Separator = 5.0 vel. heads 


Pipe friction = 1.3 vel. heads 


2 elbows = 1.2 vel. heads 
Contraction = (0.5 vel. heads 
Total = 9.5 vel. heads 
Therefore, N. = 9.5 


6. Evaporator pressure 

From the design graph for k = 1.0, and using N = 9.5 
and G/G, = 0.804, the value for P,/P,, is found to be 
0.341. Therefore, the equipment pressure will be found 
to be approximately 15 mm Hg., ie. P, = 5.0/0.341 = 
14.65 mm Hg. 


Conclusions 

1. The proposed design graphs are equivalent to | aapple’s 
as they are based upon the same assumptions. 

2. The graphs may be used for the solution of problems 
involving the flow of compressible fluids in pipes. 

3. Use of the graphs will eliminate a trial and error 
solution, previously necessary, in the calculation of the 
upstream pressure when the pipe flow rate and down- 
stream pressure are known. 


4. The graphs are especially useful in the design or 
evaluation of a complex venting system. 
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Nomenclature 


c = acoustic velocity ft./sec. 

C, = specific heat at constant volume B.t.u./(Ib.) (°F.) 
C, = specific heat at constant pressure B.t.u./(Ib.) (°F.) 
d= = derivative of 

I) = inside diameter of pipe ft. 

( = natural logarithmic base 

F = friction energy loss ft. 

f = friction factor 

g = gravitational constant ft./sec.? 

G = mass velocity Ib./(sec.) (sq. ft.) 
J) = Mechanical equivalent of heat 778 ft. lb./B.T.U. 

k = ratio of specific heats 

la = natural logarithm 

L = length of pipe bt. 


M = Mach No. 
Mw = Molecular weight of fluid 


N= Friction term 
P = absolute pressure Ib./sq. ft. 
< = gas constant = 1546 ft.lb./°F. 
= 2783 ft.lb./°C 
R. = Reynolds number 
T = absolute temperature 
l = internal energy B.t.u./Ib. 
\ = specific volume Cu. ft./Ib. 
V = velocity ft./sec. 
= = viscosity Ib./(ft. . (sec. ) 
9 = density Ib./cu. 
Subscripts 
o = in reservoir ahead of pipe 
1 = upstream point in pipe 
2 = downstream point in pipe 
C = pertains to acoustic flow 
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APPENDIX 

Derivation of Equations 
For an adiabatic expansion, the first law of thermo- 

dynamics may be expressed: 
OS me hd se Pave oaks hie hin D 

From the Ideal gas law, 
; d (Pv) , 
JdU = Gets eh eeeeeees ‘ove 


k — 


An overall energy balance results in the equation, 


V2 
O = a( ) +3004 CRs ses <n eO 
2¢ 


Elimination of Jdu from equations (2) and (3) pro- 
vides a differential form of the equation of state. 


(k - DY V? , 
k ( 2g satiate : (4) 


Integration of equation (4) yields the equation of state. 


(k — 1) 
Pivi + 


d (Pv) 


Vi? = Pove + V2?.. 34) 


Pivi = Povy + [V»2 i) on : ‘ .(6) 


By material balance 


V; Vi 
: (8) 
V2 V2 
The acoustic velocity is related to V, and M, b 
M, VkePavs (9) 








































: Bas Se 
icGraw 


itions”’, 


‘rimo- 


(%) 


.(2) 


(3) 


pro- 


(4) 


state. 


.(5) 











Substituting Eqs. (8) and (9) into (7), 


k - 1 vi’ 
Pivi = Pave + : M.? Pov. | 1 — ee) 
2 V2 
| k -1 , wi 
= Pove - + M.? | ooh, oo mS ere (11) 
| 2 V2 J 
2 k — 1 P Vi : 
yi + M.?] 1 — ne oe Oe 
P, vi | Zz Ve J 


Similarly, the equation of state based on known up- 
stream conditions can be derived. 


P, Vi k — 1 , Vo ) 
= yi + Mit b= Oe 8 ibeiscu (13) 
P; V2 l V1 J 


Eq. (12) is the desired equation which relates pressure 
and volume ratios with the Mach number based on down- 
stream conditions. Eq. (13) was used by Lapple in 
development of the design graphs based on known up- 
stream conditions. It should be noted that if k = 1, the 
equation of state reduces to P,v, = P,v,, the equation 
for isothermal expansion. 


ll 


tO 


The other equation that is required for description of 
the adiabatic flow process is the mechanical energy balance 
equation. Elimination of JdU from Eq. (1) and (3) results 
in Bernoulli’s equation, with the restriction of horizontal 


flow. 

v2 : 

| +vdP + dF =0O..... oe (14) 
2g 
Allowance for friction is provided by the relation, 
_ fV2dL ; 
GF me —@§—. on cc ccccccvcces .(15) 
2gD 


Substitution of (15) in (14) results in 


V2 dP Vidl _o 16 
( 2 + vo bt 2D = as nets (16) 


Integration of Eq. (16) yields 


Ee | i 
n v, “aes V1 os ees (17) 


Rearrangement of Eq. (17) and using 


_fL ‘4 
= Doce an One 


and 
V V 
M = == crak eae eatat 3 (19) 
© AfkgPy 
vields the following forms of the flow equation. 
+(k -—1 


) M,? V2 ; = ; 
rey | ee . (20) 
(k + 1) M2? v1 o 
+ (k — oe]: re (“) | at () an 
(k + 1) Mi? ve ie 


2| 


2kN 2 
k+1 
2kN 2 
‘+i 
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Eq. (21) was used by Lapple for construction of his 
flow graphs. The design graphs reported here were based 
on equations (12) and (20). 

The temperature relationships shown in Figures 2 and 
3 were obtained from the equation of state (12) and the 
Ideal Gas Law. The derivation is shown below. 


The equation of state is 


P re | (k — 1) M2? 3 
woes | — = ; (12) 
P. Vii 2 | Vo 


A form of the Ideal Gas Law is 


Elimination of — from Eq. (12) and (22) provides 


T,, k, and M,. 


the desired relationship between P,/P,, T, 


t— ii ese (23) 


\ oo BS 


Eq. (23) v = plotted on Figures 2 and 3 for constant 


values of T.,/T,. The minimum temperature ratio can be 

found by ae the expression inside the square root 

sign in Eq. (23) to zero, and putting M, equal to unity. 
ri/Te -— 1 

i.e. 1 — = © (24) 


Rearrangement of terms yields 


? 


r./T, = — s (25) 


k+1 


Values of the minimum temperature are shown in 
Table 1 as a function of k. 


TABLE 1 


MINIMUM TEMPERATURE RATIO 


k T2/T, 
1.0 1.000 
1.1 0.952 
RY 0.909 
3 0.870 
1.4 0.833 
Fes 0.800 
1.6 0.769 
bat 0.741 
1.8 0.714 
* * * 











High Energy Fuels for 


Jet Propulsion Engines' 


C. W. PERRY’ 


Since World War II, the development of jet pro- 
pulsion engines has fallen into two main categories, 
i.e. air breathing devices which include the turbojet 
and ramjet engines, and the rockets which use solid 
or liquid fuels. A brief review is given of the history 
and operating principle of each type. 


The use of fuel energy for power release in jet 
propulsion engines is based on rate of change of 
momentum as the fuel is burned either with oxygen 
or other oxidizer. The basic power and thrust equa- 
tions are presented. From these equations it is shown 
that high heat of combustion, AH, is the most sig- 
nificant fuel property for air breathing engines, but 
other properties including molecular weight of the 
products of combustion and temperature of com- 
bustion become equally significant for the rockets. 
Other desirable properties for chemical fuels are 
listed. 

Candidate chemical fuels having higher energy 
release than hydrocarbons (with oxygen) are tabu- 
lated. It is shown that hydrogen, beryllium and boron 
hydrides have significant potentiality. Hydrogen is 
difficult to use. Beryllium chemicals are very toxic. 
Consequently, the boron hydrides and their alkylated 
compounds show exceptional promise. The properties 


to 


of these are listed. 
ortp War Il gave tremendous the 


7 cence of all forms of jet propulsion in the 


broad sense that this includes most or all of the engines 
using the reaction principle, those that obtain their 
propulsive thrust by the discharge of a jet of fluid. The 
military need during World War II for the increased 
speed, range and altitude possible for aircraft or missiles 


impetus 


1.e, 


of 


from these engines brings to mind the statement of the 
early German military leader, who claimed 
that war is responsible for most human progress, 
another of saying that mother 
invention. 


Clausewitz, 
This is 


way necessity is the of 
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Mechanical development of jet propulsion engines 
preceded World War II to some extent. Early mechan- 
ical designs were made primarily around hy drocarbon or 
organic liquid fuels and nitrogen based solid fuels. The 
practicality of jet propulsion engines having been es- 
tablished mechanically, based on these conv entional fuels. 
it followed that added speed, altitude or range would 
be sought from higher energy fuels, including chemicals, 
free radicals, and atomic energy. High energy chemical 
fuels have been developed which are the main subject 
of this paper. The applic ation of these fuels requires 
added development of materials for engine construction 
to handle higher combustion temperatures and new pro- 
ducts of combustion. 


Jet propulsion engines fall basically into only two 
groups. The first a has been called thermal jet 
engines in the U. S. These are generally air breathing 


devices in which a stream of compressed air is mixed 
with the products of combustion from the fuel to raise 
the mixture to a desired temperature. There are two 
main subdivisions in this group, as follows: 

a. Turbojet engines 
Ramjet engines. 
Ihe pulsejet engine is, in principle, a variation of the 
ramjet engine. 

The second group of jet propulsion engines is the 


rockets, which may be subdivided into 

Liquid propellant rocket engines 

b. Solid propellant rocket engines. 
Before chemical fuels for these engines, 
the basic construction and operation of them will be 
discussed briefly; since this will give considerable insight 


and characteristics. 


discussing 


into fuel needs 


The turbojet engine 


Figure | shows the principal elements of a turbojet 
\ Now g ve Exhoust 
le ys 
Compressor —, Turbine System 


= 


Air Iniet t 


k Engine 
bs Exhaust 
Accessories and 

Accessory Drive Gear 





Figure 1—Schematic diagram of a turbojet engine. 
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engine. The modern conception of this engine can be 
trac.d to a French patent in 1922” In Great Britain 
a more specific patent was issued to Whittle in 1930. 
Phe first British turbojet engine was developed in 1937. 
In Germany, work on this engine dates back to 1930, and 


the Germans first flew an airplane propelled by a turbojet 


engine on August 27, 1939. 


Air as the oxidizing medium, is pulled or rammed into 
the air inlet depending on the forward velocity of the 
engine. It is first partially compressed by passing through 
a diverging nozzle, or diffuser, in which kinetic energy 
is partially converted to static pressure. It 1s then in- 
creased in pressure further by an axial flow or centrifugal 
compresso! At this elevated pressure it is passed through 
combustion chambers where it is heated to around 1500°F. 


(with conventional hydrocarbon fuels). 


Combustion takcs place at substantially constant pres- 
sure. The hot air, containing less than 5% combustion 
products, passes through the turbine to extract only that 
energy needed to drive the compressor. The fluid or 
gas mixture is expanded and discharged through an ex 
haust nozzie to the atmosphere. 


The ramjet engine 


The invention of the ramjet engine is credited to M. 
Lorin in France in 1913. ‘*) “) — It is, in essence, a turbojet 
engine with the compressor and turbine absent. ‘These 
Jatter components of the turbojet engine gradually become 
less effective as Hight speed increases. 


Ihe ramjet engine consists of the following elements 


? 


as shown in Figure 2 \ir entering the diffuser is slowed 


Combustion Exhaust 
Chamber Nozzle 


Diffuser 
Section 


| Exhaust 


me : 
Entrance ; Fluid 
= «SEO 


Fuel ™ 
Nozzles 


Figure 2—Schematic diagram of a ramjet engine. 







down 1n it to subsonic velocity to convert kinetic energy 
into static pressure. Then it passes into the combustion 
chamber where its temperature ts raised by burning a fuel. 
Phe hot gas mixture in the combustion chamber then ex- 
pands through a nozzle and is discharged to the atmos 
phere. 


[he ramjet engine has the advantage of few moving 
parts, and consequently an apparent simple design. In 
practise it is difficult to design and develop, however. It 
must be placed in operation at a high speed, which re 
quires either rocket launching or dropping from an air 
plane. Its best area of operation is from the upper 
speed limit of the turbojet engine (Mach 2 roughly) to 
about Mach 4. It has a definite place in the overall 
£ Advanced de 
signs have been developed in the U.S. by the Marquardt 
\ircratt ( omipany in California 


development of jet propulsion engines. 
I pro} 


\ variation of the ramjet engine is the pulseyjet, which 
Lhe Germans 
vorking engine during World War IIL to 


powe! the V-1 buzz bombs 


is in essence an intermittent flow ramjet. 


developed a 
Through a valve system, 
the burning of fuel takes place in the combustion 
chambers from 40 to 250 times per second. It has the 


advantage ot developing thrust at zero speed, and is well 
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adapted to subsonic unmanned aircraft. Workable pulse. 
jets have been developed by American and British com- 
panies. 
Rocket engines 

As previously stated, rocket engines are classified into 
two types depending on the propellants. The fuel and 
oxidizer may both be liquids, in which case the engine 
is termed a bipropellant liquid one. In the second case, 
the fuel and oxidizer are blended together (either chem- 
ically or physically) as a solid mass which will burn uni- 
formly along a given axis to generate a constant pressure. 
This is termed the solid propellant rocket engine. 

In its simplest terms, a liquid bipropellant rocket en- 

gine consists of the following elements: 

1. A combustion chamber where fuel and oxidizer 

react, creating a high pressure and temperature gas 

mixture. 
2. A diverging nozzle for expanding the combustion 
gases and providing thrust. 

3. A propellant storage system consisting of separate 
tanks for fuel and oxidizer. 

4. A propellant feed system consisting of turbine 
driven pumps or a gas pressuring system, and pro- 
pellant meters and controls. 

A common configuration of these elements is given in 
Figure 3, 


Turbine Propellant Cumbentien 
Gas Nozzles Soaks 
Generator samy 
\ Jet Nozzle 
Payload — | \ J 


rl l{ 


Turbine 





Guidance » seed 
Fuel Propellant \ “Regenerative 
Oxidant Injection Cooling 

Pumps 


Figure 3—Schematic diagram of a liquid propellant 
rocket engine. 





Che solid propellant rocket engine is much simpler 
than the bi-propellant liquid one in that the fuel and 
oxidizer are either a physical solid mixture or a mole- 
cularly combined compound placed in a storage tank that 
becomes the combustion chamber as the fuel and oxidizer 
react. The hot gas mixture at high temperature and 
pressure is allowed to expand through a nozzle to provide 





thrust. A common configuration is shown in Figure 4. 
< Payload \ Propellants < Jet Nozzle 
7 
( | 
had 
/ 
Guidance Controls “_ Combustion 


FKigure 4—Schematic diagram of a solid propellant 
rocket engine. 
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Rocket engines provide the only practical means now 
available for space travel beyond the earth’s s atmosphere. 
Consequently they are now receiving extraordinary at- 
tention. It is commonly accepted that the liquid bi- 
propellant rocket engine has the best chance for long 
range travel. The rate of fuel usage can also be con- 
trolled. The solid propellant rocket engines are simple, 
less expensive, and quite reliable. 

Since rocket engines of either type must carry their 
own oxidizer, their rate of fuel consumption is very large. 
A 10,000 pound thrust rocket engine, for ex ample, will 
use 45 pounds of propellants per second. The time of 
a flights ~~ all present- ped rockets, is quite short. 
The German V-2 missile had a powered flight time of 
67 seconds, as an example. 


Basic thrust and power equations for thermal jets 

All jet propulsion engines operate on Newton’s second 
and third laws of motion. ‘2) For the thermal jet, or air 
breathing engines, the weight (pounds) of material flow- 
ing through them per unit of time (seconds) is: 


(lb./sec.) = Gair + G fuel... 1 
It follows that the force or thrust, F, is 


Ga + Gi G.V 
F (Ib.) = w - + (Pe + Po) Ac....+..(2) 
gS g 
In this equation, q is the gravitational constant, w is 
the relative jet velocity (relative to the exhaust nozzle) 
(ft./sec), V is the flight velocity A, is the entrance 
cross section, P, and P, are the entrance and exit static 
me G.V 

pressures. The term is the “ram drag” which 
applies only to the air entering, and (p. - p,) Ae is the 
“pressure thrust”. This latter term is usually small 
compared to the momentum thrust. Since most thermal 
jet engines (turbojets and ramjets) are continuous flow 
devices, this force is continuously applied to the plane 

or missile. 


Eq. (2) may be simplified by neglecting the pressure 


thrust to 
; 1 GV 5 
F = 1 - (9) 
\ g 


where v = speed ratio, V is flight velocity and 


w is the relative jet velocity. 
To convert force F in pounds to power produced 
from a thermal jet engine the heat of combustion of 


the fuel must be considered. Consequently 


\" 
G {| J4 H+ : ‘ (4) 


where P is power input, Gy is fuel feed rate (lb./sec.), J is 
the mechanical equivalent of heat, H is the heat of com- 
bustion, and V is flight velocity. 


P in (ft. Ibs. /sec 


Another way of examining the required characteristics 
for fuels for thermal j t engines is to examine the basic 
range equation, sick) is as follows: 


I N (AH) (L/D) log (W W,) (5) 


where R is the range of an aircraft, N, is the propulsion 
eficiency. AH is the heat of combustion of the fuel, 
(L/D) is the ratio of lift to drag, and W,/W, is the 
ratio of the weight of fuel to gross weight of the aircraft. 

This indicates that the range or length of a flight is 
directly proportional to AH. Pertinent factors in high 
energy liquid fuels for thermal jet engines are: 

l. High heat of combustion 

2. High density 
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3. Low volatility 

4. Chemical stability 

5. Low freezing point 

6. High boiling point 

7. Low viscosity 

8. Availability © 

9. Cost of production 

Needless to say the thermal jet engines have been 
developed around liquid fuels which may be throttled 
during flight and the fuel replaced for repeat flights. 


Basic thrust and power equations for rockets 

In the case of a rocket engine, the fuel is oxidized 
with a self-contained chemical. Therefore, “total= © 
oxidizer + © fuel, and 


Giees 
Fins -( )s + (Pe — Po) Ac. S64 QGee 
g 


This may be simplified to 


g 


where V, is the Effective Exhaust Velocity. Com- 
bining the above equations gives 


Qo 
S 


V; = w + & (Ps — Peles coc erin (8) 


The use of Effective Exhaust Velocity, (a fictitious 
velocity), is convenient because in testing rocket en- 
gines, F and G are readily measured. The ratio F/G 
is called the practical “specific thrust” or “specific im- 
pulse”. Thus the practical specific impulse is 

IK ot 
foi 05 ON Er icla i SR (9) 
G - G; xg 
If F is in pounds, and G is in pounds per second, then 


l,, isin seconds. 
Where the nozzle of the rocket motor expands the 
gases completely, so that P, = P,, then V,; = V exhaust. 
For rocket engine fuels new factors come into play. 
From the basic equations just developed, it is apparent 
that AH for the fueleoxidizer combination is of direct 
importance. For liquid rocket fuels, all of the factors 
mentioned for thermal jets are important. For either 
liquid or solid rocket fuels, since the effective exhaust 
velocity, V, of the combustion products (from fuel and 
oxidizer) is of critical importance, the average molecular 
weight of these gases becomes significant. Specific im- 
pulse may be ¢ -alculated on a theoretical basis by thermo- 
chemical methods. In this case the theoretical impulse 


is 


I.» (theoretical) = K k te 7, -. (10) 


Vk +1 \M 


where K is a proportionality constant, k is the ratio of 
specific heats, t. is the temperature of combustion, and M 
is the average molecular weight of the combustion pro- 
ducts. Z, is an expansion factor. In examining candi- 


date rocket fuels, the factor is of prime importance. 


t 
Vn 
\greement between practical and theoretical I,, 
from .93 to .98 for well designed engines. 


ranges 


For a rocket motor, where G, and G, are added to 


give G then 
' - 
Pin = (Go + Ge) JAH+> atdceeadasOnee 
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This indicates the fundamental importance of a high 
heat of combustion in the fuel, or fuel-oxidizer combina- 
tion. 

Rocket motors generally operate between 3000°F. and 


8000°F., and 200 to 1500 psia. Consequently, in cal- 
culating theoretical I,,, it is gener rally assumed that 
1. The isentropic expansion process is complete from 
combustion pressure (p,) to external pressure (p,) 
2. The equilibrium composition of gases in the com- 
bustion chamber is unaltered during the nozzle ex- 
pansion (frozen equilibrium). 


For liquid rocket propellants the following factors are 
most important: 

1. High heat of combustion (for the fuel-oxidizer 

system ) 

2. High reaction rate 

3 High propellant density 

4,6 hemical stability 

5. High boiling point 

6. Low freezing point 

7. Availability 

8. Cost of production. 

With small variations the same factors apply to solid 
rocket propellants. 


Candidate high energy fuels for jet propulsion engines 


Consider now the candidate high energy fuels for 
jet propulsion engines. For comparison, our examination 
of present fuels is necessary. 


Most of the thermal jet engines today are fueled with 
petroleum fractions overlapping the gasoline and kerosene 
volatility ranges. As such, they have a heat of com- 
bustion of about 19,000 B.t.u. per pound when burned 
with air. Special hydrocarbons having a high content of 
condensed napthenic rings may be found w hich can give 
about 5% more B.t.u. per unit volume of fuel. These 
have not been developed except experimentally, and have 
shown more success for the purpose of chemical stability 
than added energy. The most common liquid rocket 
fuels are hy drocarbons or ethyl alcohol burned with liquid 
oxygen, nitric acid or hy drogen peroxide as the oxidizer. 


Solid rocket propellants are of two types, both of 
which must combine fuel and oxidizer, either physically 
or chemically. The first, or composite type, to date has 
usually consisted of organic poly mers or rubbers as the 
fuels, mixed physically with nitrates or perchlorates as 
the oxidizers. Common oxidizers are ammonium nitrate 
and ammonium perchlorate. Additives and inhibitors are 
used to control burning. 5) The second type contains 
chemical compounds such as nitrocellulose or nitrogly- 
cerine which will support combustion by themselves. 
These two chemicals are often combined as colloids, and 
are referred to as double base solid propellants. 


Figure 5 gives the heats of combustion of several 
present or candidate fuels for jet propulsion engines. JP-4 
is the conventional hydrocarbon jet fuel used for turbojet 
engines and rockets. (6) Note the relativ ely low heats for 
ethy! alcohol and hydrazine, which are, nevertheless, good 
rocket fuels. This illustrates the point that AH does 
not have the same unique importance in rocket fuels as 
it does for the thermal jets, and low molecular weight of 
combustion products is of direct importance. 

Hydrogen gas with its large heat of combustion of 
51,570 B.t.u. /Ib. with oxygen would appear to be of ex- 
ceptional promise as a fuel. For the practical reasons of 
low liquid density, refrigeration problems (from its low 
boiling point), and difficulty i in storage and handling, hy- 
drogen i is not in use as a fuel for jet propulsion engines. 
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Beryllium metal at 29,140 B.t.u./lb, and its hydride at 
33,000 B.t.u. /\b. with oxygen also appears to have excep- 
tional promise. Again for practical reasons beryllium 
compounds are limited by extreme toxicity, limited avail- 
ability, and difficult metallurgy. 


Boron fuels 


The high energy fuels considered to be of exceptional 
promise currently in the United States for both thermal 
jet and rocket engines are the boron hydrides and their 
derivatives. An examination of the heat of combustion 
of all of the common metals and some of their hydrides 
indicates that, next to hydrogen and beryllium, boron is 
most attractive. This is substantiated by extensive re- 
search and development conducted during the last several 
years by Olin Mathieson Chemical Corporation and others 
under U. S. Navy and Air Force sponsorship. This work 
has resulted in placing two semi-production plants in 
operation and a tonnage plant being built at Niagara Falls, 
N.Y. by Olin Mathieson, the latter for the Air Force. 


These boron fuels are being developed on a sound 
basis of adequate ore reserves within the continental 
United States, and proven methods for converting borax 
(Na,B,O,) or other boron ores to usable commercial 


compounds such as B,O,, H,BO, and B,C. 


Boron chemistry 


From the basic boron chemicals mentioned above, 
more active intermediates such BCI, and NaBH, can 
be prepared. These two chemicals are now in commer- 
cial production in the U.S. From these the simplest boron 
hydride, diborane (B,H,) can be prepared. One method 
is as follows (7): 


2NaBH, + 2HCl — B2H, + 2NaCl + 2H: 


Other acids may also be used with sodium borohydride to 
give diborane and the corresponding salt. 


A second synthesis of diborane, which has been used 
in the U. S. is as follows: (8) 


6LiH + 2BF; (C2H;).0 -}B2H, + 6LIF + 2 (C2H;).0 


The boron trifluoride etherate is available commercial- 
ly in the U.S. Other metal hydrides and other boron 
halides may be used in analogous manner. 


Diborane may be converted to higher boron hydrides 
with lower hy drogen contents by pyrolytic reactions as 
follows: (9 

B,H, <— 2BH, 

BH, + |B, H, — Intermediates 

B, H, a Intermediates — B,HI,, + 2H, 

B. H, — B.H, + H. 


11 


B.H,, a B, H., + Higher Hydrides 
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The boron hydrides range from the gaseous B:He 
through liquids to solid decaborane, and higher solids. 
The properties of three of these hydrides are as follows: 


PHYSICAL PROPERTIES OF BORANES 


Mol. M.P. B.P. Sp. 

Wt. “Ps a Gr. 
Diborane(B:H,) | 27.7 | -265 | —134 | 0.43 
Pentaborane (B;Hg) 63.2 | — 52 140 0.61 
Decaborane (BioH1s) 122.3 211 415 0.94 (s) 


These boron hydrides (excepting decaborane ) are 
somewhat difficult to handle because of their toxicity and 
pyrophoric nature. They may be modified by adding 
alkyl groups so that they become stable liquids however, 
with low toxicity and other satisfactory physical proper- 
ties. These liquid fuels are little more difficult to store 
and handle than leaded gasoline. The heats of combus- 
tion of the boron hydrides, and some of their alkyl de- 
rivatives are given in Figure 6. 
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Figure 6—Heats of combustion, B. hydrides. 


Application of the boron fuels 

Stable liquid and solid boron fuels having high heats 
of combustion and other desirable physical properties are 
of interest for use in turbojets, and their afterburners, 
ramjets, and rockets. The small quantities of these boron 
fuels now becoming available substantiate the promise 
that they hold. The exhaust products from air breathing 
engines containing B,O, and water vapor of low toxicity. 

For the thermal jet or air breathing engines the writer 
considers that the modified boron hydride liquid fuels 
offer the best promise when considered from the overall 
point of view of availability of raw materials, manufac- 
turing technology, physical and chemical properties, and 
production cost. Since the higher boron hydrides are 
liquids and solids, they can also be considered as candidate 
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fuels for bipropellant liquid and composite solid rocket 
engines. 


Other fuels 


In the case of rocket engines the fuel and oxidizer 
must be considered in combination rather than separately. 
Ratios of oxidizer to fuel may be either equal to or greater 
than the stoichiometry requires, and this of course affects 
the specific impulse. It has been shown that high heat of 
combustion for a rocket fuel is of interest primarily to 
obtain a high combustion temperature, and that low mole- 
cular weight of combustion product is also of prime im- 
portance. Figure 7 gives theoretical, or calculated speci- 
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Figure 7—Theoretical specific impulses of rocket 
Fuels and oxidizers. 


fic impulses for present and future fuel-oxidizer combina- 
tions, as well as free radicals. '°) Theoretically, at least, 
free radicals offer possibilities beyond hydrogen ‘and ozone 
for high specific impulses. This is far from a reality at 
the present time, however. 

Note that the liquid bipropellant combinations cur- 
rently in use for rockets, ethyl alcohol-red fuming nitric 
acid, kerosene (JP-4) - 99% hydrogen peroxide, etc. all 
give I,,, below 250 seconds. The better oxidizers, ozone 
and fluorine, offer attractive increases in I,,, but are hard 
to handle. Note that diborane and fluorine given an I,, 
of 291. All of these combinations are figured at 300 psia. 


Conclusion 

It is safe to say that jet propulsion engines will pre- 
dominate as power sources for aircraft and missiles in the 
immediate future. They will to an increasing extent be 
propelled by high energy chemical fuels to give greater 
speed, altitude, load and range. These chemical fuels will 
probably in the future, be supplemented by free radical 
and atomic energy sources. The chemical fuels will have 
a permanent place in the propulsion of military and com- 
mercial aircraft and missiles, however. 
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Mass Transfer in a Bubble 


Column’ 


A. I. JOHNSON? and C. W. BOWMAN? 


The effects of gas and liquid flow rates and of 
seal height on the absorption of carbon dioxide in 
water have been studied. The gas bubbles were 
formed from 1/16 in. orifices, with both single and 
multiple hole arrangements. 


Bubble frequencies were measured using a high 
speed motion picture camera and are reported graph- 
ically as a function of gas flow rate. At high flow 
rates, considerable breaking up of the bubbles was 
observed. 


Because of the difficulty of estimating interfacial 
areas, the transfer data were reported using the 
product of the transfer coefficient and the interfacial 
area, K,a. For single orifices K,a was found to vary 
directly with the gas rate, was only slightly dependent 
on liquid rate and was almost independent of seal 
height after due allowance for end effects. With 
multiple hole plates having °% in. spacing, there was 
no appreciable change in the transfer rate per hole 
until the hole density exceeded 1.5 per square inch, 
after which the rate decreased slightly. 

End effect contributions to transfer are also 
reported. 


Hr study of the rate of transfer of material between 
Lae fluid phases is of high theoretical and practical 
significance. By analogy to the equation for the rate of 
diffusion in gases, it is usual to define a rate equation for 
diffusion in liquids of the form: 


N=K,, (x,—x) 


Several theories have been advanced to relate K,, to 
the diffusivity in the liquid phase. The steady state 
stagnant film theory “) predicts K;, to vary directly with 
the diffusivity. However theories om on unsteady state 
transfer indicate that K,, varies as the square root ‘of the 
diffusivity 2%). In ectual practice, K;, is found to vary 
as 1.5 to 1D%-6 (2. 5) which lends support to the unsteady 
state transfer mechanism. 


Bubble type vapor-liquid contactors have frequently 
been used to study the relationship between K,, and 
liquid diffusivity, viscosity, density, surface tension, tem- 
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perature, and the flow properties of the system. How- 
ever the evaluation of interfacial area in such apparatus 
presents considerable difficulty. Many investigators have 
assumed that the rising bubbles are spheres and have thus 
been able to calculate the surface area available for trans- 
fer. This is a questionable approach in the low gas 
throughput (small bubble) region and is completely un- 
tenable at higher rates. A statistical analysis of the bub- 
ble size and shape distribution seems necessary. 

This research was a study of the effects of the flow 
variables on the group “K,a” in a bubble column at high 
gas throughputs. The bubbles were generated from 
1/16” diameter holes into a water seal varying from one 
to six inches. The number of transfer units was cal- 
culated based on complete vertical and horizontal mixing 
in the continuous phase (i.e. — uniform concentration). 
K,a was evaluated from the equation: 


CLé 2S 
(N.T.U.)ou = Rise 
4L 


Theoretical principles 


(a) Transfer unit concept: 


For a gas-liquid countercurrent contactor, the follow- 
ing equation can be written: 


bes Kaa (t= 2)occsccik cess 


dx K,aV 
ae eee er TEER E TCT CT Tt (2) 
Xe — X L 


The left hand side is defined as the number of over- 
all transfer units based on the liquid driving force. For 
the cylindrical column used in this study, the above 
equation can be expressed as follows: 


Xt 


Pl dx K,amrd2s 
(N.T.U.)or = - sarees wean 
Xe —X 4L 


X2 


(N.T.U.)o;, can be determined experimentally—thus per- 
mitting evaluation of K,a. 


(b) Measurement of the number of transfer units 


In the present investigation, the liquid (water) in the 
column was assumed to be of uniform concentration. (This 
was indicated by observing the dispersion of dye solu- 
tions injected at the top of the column.) For this case, 
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the expression for (N.T.U.) 5; reduces to where E is 


E 
1-E 
the efficiency (fractional approach to equilibrium) de- 


ig . Pe ed ee : 
fined by the relationship — 2. Since x, was zero, this 


xX, 

. a x 
simplifies to —* 
Xe 

For the case where a concentration gradient exists in 

the column, (N.T.U.)o,, is given by-In(1-E). This latter 

expression was not used in this paper except for com- 

parison purposes. 


(c) End effect correction 


Transfer occurring during the formation of a bubble 
and also during its collapse at the top of the column is 
not by the same mechanism which occurs during the 
rise period. It is convenient, therefore, to treat this 
anomalous transfer separately under the term “end 
effects”. 

The usual method of determining the magnitude of 
this abnormal terminal transfer is to plot the number of 
transfer units for the over-all process against seal height 
and to extrapolate to zero seal height. The intercept gives 
an approximation of the number of transfer units due to 
end effects. 

The number of transfer units in the steady rise part of 
the column for any seal is calculated by subtraction of 
this intercept value from the total number of transfer 
units. 


By treating the number of transfer units for the end 
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Figure 2 


effect using the ‘ ‘completely mixed” relationship, the cor- 
responding efficiency E,, is obtained. 


(d) Literature review 


Several workers studying mass transfer from single 
bubbles of a fairly insoluble pure gas rising in water have 
obtained K,, values from 100 to 200 cm/hr. (7), The 
bubbles were assumed to be spheres. 


Two investigators ‘> 6) reported lower values of K,, 
for bubbles below a certain size. Garner and Hammer- 
ton) found that for bubbles less than 0.2 cm. in diameter, 
the transfer coefficient dropped to about 25 em/hr. 
which is the value predicted from the Fréssling equa- 
tion’) for diffusion from rigid spheres. The departure 
from rigid sphere behavior occurring at a critical diameter 
was atttributed to the onset of internal circulation in the 
bubble. 


Two studies (1% 11) have been carried out on mass 
transfer from bubble swarms at very high gas rates. The 
results were expressed as a “K,,a” product or its equivalent. 


Many investigations have been undertaken on the 
mechanics of bubble behavior in liquids. The work has 
been concerned chiefly with (a) the measurement of 
bubble volume (or frequency ) as a function of the gas 
rate through an orifice and (b) the velocity of rise as a 
function of bubble diameter. An excellent theoretical 
and experimental treatise has been given by Siemes in a 
series of papers(12 13, 14), 

Quite recently some workers '* 15) have pioneered an 
attack on the problem of analysing bubble behavior at 
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high throughputs where a statistical approach is neces- 
sary. Bubbles at these high gas rates can no longer be 
considered as spheres. It has been found convenient to 
define three bubble categories—spheres, ellipsoids, and 
flat bottomed domes. 

Unfortunately, studies on mass transfer have not kept 
yace with this rapid development in bubble mechanics. 
The ideas of Siemes, particularly in his paper on the 
“Dispersion of Gases in Liquids by Means of Nozzles 
at High Throughputs”'*), should now be applied to 
systems involving mass transfer. 


Experimental details 


See Figure | for a diagrammatic sketch of the apparatus: 
(the detail of the lower portion of the bubble column is 
shown in Figure 2). 

Briefly, carbon dioxide from a supply reservoir (14) 
was bubbled continuously through a water seal (5) until 
a definite volume had been absorbed, as indicated by the 
change in level in water reservoir (19). The water flow 
to and from the column was maintained at a constant 
value by the automatic control valve (3) and the needle 
valve (7). 

To calculate the efficiency of the column, the follow- 
ing information was required for each run: 


1. Volume change in reservoir (19) 


2. Leakage rate of gas from the system 
3. Time of run 
4. Temper ature of water in bubble column (5) 


(controlled to 25+ 0.1°C.) 
Temperature of gas in reservoir (14) 
(range — 23-26°C.) 
6. Water flow rate 
Average system pressure. 


wu 


In addition the following information was recorded: 
1. Seal height 
2. Gas rate 
3. Temperature of gas in surge flask (10) 
4. Water level in reservoir (14) 


The experimental investigation was divided into four 
groups as follows: 


Group I 


Exploratory runs using raw water from the city mains. 


Group II 


Effect of gas and liquid flow rates on transfer rate 
using distilled water. (seal height constant at six 
inches). See Figures 3 to 10. 


Group Ill 


Effect of seal height on transfer rate using distilled 
water. See Figures 11 to 16. 


Group IV 


Multiple hole orifice studies. See Figures 17, 18. 
All results are tabled in Reference 16. 
Discussion of results 


(a) Frequency measurements 


Before proceeding with an analysis of the mass transfer 
data, the results of a high speed motion picture study on 





The Canadian Journal of Chemical Engineering, December, 1958 


bubble frequency (with 0.159 cm. orifice diameter) will 
be given. 

At low gas rates (below 5 cc./sec.), the bubbles formed 
and rose singly without physically contacting one another. 
At slightly higher throughputs, a bubble did not have 
time to leave the vicinity of the orifice before the next 
bubble had started to grow. At some critical flow rate, 
the forming bubble w ‘ould actually penetrate its prede- 
cessor with its interface remaining intact for a brief 
moment. These two bubbles would then merge and 
rise as a single unit. At still higher gas rates, three and 
sometimes four bubbles would merge in this fashion be- 
fore a break occurred in the chain. Above 50 cc./sec., 
an almost continuous jet emanated from the orifice with 
the bell-like top breaking away at fairly uniform time 
intervals. 


Thus up to 5 cc./sec., the frequency as measured at 
the orifice was identical to that measured at any other 
point in the column. However, above this rate, the fre- 
quency as observed even a few centimeters above the 
orifice suddenly decreased to less than half its previous 
value due to the coalescence of successive bubbles. The 
frequency of these bubble “clusters” remained fairly 
constant over the rest of the throughput range. Figure 
19 illustrates this phenomenon. 

As far as mass transfer studies are concerned, it is the 
behavior of the bubble cluster during the rise period that 
is of main interest. Above approximately 20 cc./sec., 
this large bubble mass assumed a mushroom-like shape 
which tended to disintegrate about two inches above the 
orifice This shattering effect became more pronounced 
as the gas rate was increased with smaller and smaller 
bubbles resulting. Measurement of interfacial area would 
obviously have to be carried out by a statistical approach 
similar to the method of Siemes and Gunther(!*), 


(b) Mass transfer results 


Initially, raw water from the city mains was used 
in this study but it was soon observed that the repro- 
ducibility of column efficiency measurements was very 
poor. The results obtained, therefore, will not be dis- 
cussed further. 


Group II 


In order to rule out the possibility that changes in 
raw water supply were responsible for the above noted 
variability, it was decided to repeat the above runs using 
distilled water. A marked improvement in consistency 
immediately resulted. Therefore distilled water was used 
in all subsequent test work. 


Figure 3 is a graph of measured column efficiency 
(Ey) versus gas rate with liquid rate as parameter. 
“Smoothed” values of E;, were taken from this graph 
and used to calculate (N.T.U.); (Figure +). The same 
data was plotted in Figure 5 to show the effect of liquid 
rate. 


It can be seen from these graphs that (N.T.U.),7 is 
directly proportional to gas rate (above 16 cc/sec.) and 
inversely proportional to liquid rate. This proportional- 


ity between (N.T.U.), and ¢ would be expected con- 


sidering Eq (3). The effect of gas rate is obviously the 
reflection of increasing interfacial area. 

However, as discussed previously, all such correla- 
tions should be based solely on the transfer occurring 
during the steady ‘rise’ portion of the column. There- 
fore the above data was replotted (Figures 5, 6) after 
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Figure 3—Column efficiency (with no end effect correc- Figure 4—Number of overall transfer units (with no end 
tions) versus gas flow rate for carbon dioxide bubbles effect corrections) versus gas flow rate. 
formed from a 1/16 inch orifice. 


allowing for the transfer due to end effects as deter- responding to the constant frequency region) and there- 
mined from Group III tests. The end effect corrections fore the data in Figure 6 cannot be extrapolated to low 
were only applicable to gas rates above 27 cc/sec. (cor- gas rates. 
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Figure 5—Number of overall transfer units (with and Figure 6—Number of overall transfer units (with end 
without end effect corrections) versus liquid rate. effect corrections) versus gas flow rate. 
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It can be seen that the general form of the dependence 
of N.T.U. on gas and liquid rates is not appreciably dis- 
turbed by the allowance for end effects. 

For comparison purposes, the difference between 
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Figure 9—Mass transfer coefficient (K,a—with and with- 
out end effect corrections) versus gas rate. 
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Figure 10—Mass transfer coefficient (K, a—with and with- 
out end effect corrections) versus liquid rate. 
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Figure 11—Number of overall transfer units (without end 
effect corrections) versus seal height at a carbon dioxide 
rate of 50.8 cc/sec. 
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Figure 13—Number of overall transfer units (without end 
effect corrections) versus seal height at a carbon dioxide 
rate of 26.6 cc/sec. 
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Figure 12—Number of overall transfer units (without end 
effect corrections) versus seal height at a carbon dioxide 
rate of 89.3 cc/sec. 
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Figure 14—FEfficiency contributed by end effects versus gas 
flow rate. 
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Figure 15—Mass transfer coefficient (K,a—without end 
effect corrections) versus seal height. 
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Figure 18—Number of overall transfer units per hole (for 
both single and multiple hole orifice plates) versus gas 
flow rate. 


Figure 17—Number of overall transfer units per plate (for 
both single and multiple hole orifice plates) versus gas 
flow rate. 
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FREQUENCY 
MEASUREMENTS 
BY W.SIEMES (12) 


FREQUENCY MEASUREMENTS 


Orifice diameter 0.159 cm. 


o 
9 


Orifice diameter 0.150 cm. 


S 
9 


Frequency 


O—Measured at orifice 
O—Measured 2" above orifice 


Gas rate- cc/sec 


Frequency — No. bubbles /sec. 


40 60 
Gas rate cc/sec 


Figure 19—Bubble frequency (for 1/16 inch—0.159 em— 
orifice) versus gas flow rate. 


By plotting the above data in the form of K,a versus 
gas and liquid rates (Figures 9 and 10), it was found that 
K,a does depend slightly on liquid rate for high gas 
throughputs. (Kj, values could have equally as well 
I , : 
I relationship 
was best suited for the N.T.U. graphs and the liquid rate 
was therefore carried through in this form.) Thus NCE. 
under these conditions is slightly less than directly pro- 


been plotted versus L. However, the 


portional to ; (the upper curve in Figure 5 is too steep 


to show this). This may be due to experimental error 
but since this phenomenon was observed in subsequent 
test work, it is felt that the effect was a real one. 


A reason for this increase in K,a with increasing 
liquid rate is not hard to find. The motion picture studies 
revealed that a more complete bubble break-up occurred 
at higher liquid rates. This resulted in smaller bubbles and 
hence larger values of ‘a’ 


Group III 


This series of tests was a study of the functional 
dependence of (N.T.U.), on seal height. (Figures 11, 12, 
13, 14). Three different gas rates and three different 
liquid rates were combined in all possible combinations— 
making a total of nine sets of experiments. 

As can be seen ee Figures 11, 12 and 13, the rela- 
tion between (N.T.U.), and seal height i is ebeceed by the 
flow properties of the system. Genet rally at high gas and 
medium liquid flow rates, (N.T.U.)> v varied linearly with 
seal height. At low gas rates, curves convex downward 
were obtained. 

By wrt the data back to zero height, the 
intercept (N.T.U.),- was obtained. This permitted evalu- 


ation of FE, using the 


E 
E relationship. 


E, is plotted versus gas rate with liquid rate as para- 
meter in Figure 14. Surprisingly, E, was found to be 
almost independent of gas rate. Note that in cao 13, 
for low gas and high ‘liquid rates (bottom curve), the 
total column transfer was due to end effects up to a : seal 
height of about 3 in. 


However it must be mentioned that the end effect 
was probably not independent of seal height for these 
tests. At low seal heights, the swirling motion of the 
liquid around the orifice was quite violent and the 
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distortion of the bubble during formation much greater 
than was observed with higher seals. Thus although 
Figures 11, 12, 13 do show abnormally high transfer at 
low seals, application of the “apparent” end effect as a 
correction for high seal studies may not be completely 
valid. 

For this investigation, however, the observed depen- 
dence of (N.T.U.), on gas and liquid rate was used to 
calculate (N.T.U.),,; values for the tests in Group II. 

Figures 15 and 16 show both (K,a), and (K,a)y 
values as a function of seal height. (K;,a)y, was found to 
be almost independent of height which tends to support 
the above treatment (K,a), values rise abruptly at low 
seal heights, reflecting the abnormal transfer behavior at 
bubble formation and collapse. 


It can be seen that the slopes of curves 1, 2, and 3 (high 
gas rate) in Figure 16 are in order of decreasing liquid 
rate. Curve 1 (the highest liquid rate) has a slight positive 
slope; curve 2 (medium liquid rate) is horizontal; curve 
3 (lowest liquid rate) slopes downwards. This is again 
evidence that high liquid rates tend to cause a more 
complete shattering of the bubble dome and hence larger 
interfacial areas. 


Group IV 


The multiple orifice studies were carried out on 1, 2, 
3, 5 and 7 hole plates. Refer to Figure 2 for a layout of 
the distribution of holes. 

The first tests were carried out with a constant liquid 
rate per hole so that (N.T.U.), values could be compared 
on a “per plate” basis. Thus it would be expected that 
(N.T.U.), values would be the same for all plates provided 
the bubbles from adjacent holes rose without interference. 

It can be seen from Figure 17 that the transfer rate 
increased considerably w ith an increasing number of 
holes would magnify or diminish this effect. For al- 
crease may have been due to the increasing liquid rate 
rather than actual interference among the bubbles. 


For this reason, a second series of tests was under- 
taken with a constant liquid rate per plate. The (N.T.U.)+ 
values were therefore expressed on a “per hole” basis 
to permit easy comparison of the results. Figure 18 shows 
that increasing the number of holes decreases the transfer 
rate but only to a very slight extent. (Figure 18 was 
felt to be more indicative of “bubble interference” phe- 
nomena than Figure 17). 


It is difficult to say if decreasing the spacing between 
holes would magnify or diminish this effect. For al- 
though coalescence at the orifice would be increased, 
the off-setting effect of more complete bubble break-up 
higher in the column would probably result. 


Conclusions 


As the flow of gas through a circular 1/16 in. diameter 
(0.159 cm.) orifice was increased, the bubble frequency 
rose to a maximum value (about 32 bubbles per second) 
and then suddenly decreased to a fraction (4 to 4) of this 
figure due to coalescence of successive bubbles. The 
frequency of the resulting “bubble clusters” remained 
constant over a wide range of flow rates (5 to over 
100 cc./sec.). Above 20 cc./sec., the bubbles assumed a 
mushroom shape and disintegrated after rising a few 
inches. This shattering effect became more pronounced 
as the gas rate was increased with a shift in the bubble 
size distribution curve towards smaller bubbles. 


For single orifices, K;,a varied directly with gas rate, 
was only slightly dependent on liquid rate and was al- 
most independent of seal height if allowances were made 
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for end effects. For high gas throughputs, K,a was 
found to increase slightly with increasing liquid rate 
due to a more complete break-up of the bubble “mush- 
room . 

The efficiency (fractional approach to equilibrium ) 
contributed by end effects decreased with i increasing liquid 
rate but was almost independent of gas rate. The latter 
suggests that the water in the immediate vicinity of the 
orifice approached saturation at comparatively low gas 
rates. 

The interaction among bubbles on a multiple orifice 
plate (with & in. hole spacing) was not appreciable until 
the hole density exceeded 1.5 per square inch. Above 
this figure, the ‘transfer rate per hole decreased slightly. 


With the shallow seals used in this study, the mixing 
action of the gas and liquid flows was sufficient to result 
in a uniform concentration in the continuous phase. This 
was indicated by observation of the dispersion pattern 
of dye solutions injected into the column and was later 
verified by the linear relationship obtained between 
(N.T.U.)o;, (when calculated for complete mixing of the 


. ] 
continuous phase) and L 


4 


Nomenclature 

a = interfacial area per unit volume. 

d, = column diameter — inches. 

Di = liquid diffusivity — ft.?/hr. 

Ex = overall column efficiency. 

Em = column efficiency after allowing for transfer due 
to end effects. 

Er = efficiency due to end effects. 
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G = gas flow rate — cc./sec. 
= are ae l 
Ki = overall liquid film transfer coefficient. = ft./hr. 
L = liquid flow rate — ft.?/hr. 
a = mass transfer rate per unit area. 
(N.T.U.)or = number of overall transfer units based on liquid 
a driving force. 
(N.1 -U.)r = (N.T.U.)ot based on Er. 
(N.T.U.)m = (N.T.U.)ot based on Em. 
(N.T.U.)e = (N.T.U.)or based on Er. 
S = seal height — inches. — 
\ = volume of column — ft.* 
x = concentration of CO» in water — lb./ft.* 
X1 = concentration of COs in water at bottom of column 
Ib./ft.9 
X2 = concentration of CO» in water at top of column — 
Ib./ft.3 | 
Xe = concentration of CO, in saturated water — lb/ft. 
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Recent Developments 


in the Manufacture of Chlorine Dioxide 


W. HOWARD RAPSON- 


Fundamental studies on the kinetics of formation 
of chlorine dioxide have led to an hypothesis con- 
cerning the mechanism of the reactions involved. This 
has lead to improvements in the manufacture of this 
important bleaching and oxidizing agent. Older pro- 
cesses have been made more efficient and some new 
ones have been developed. 

A new process, patented by the author and 
assigned to Electric Reduction Co. of Canada Limited 
in Canada, and to Hooker Chemical Corporation in 
the United States, employing sodium chlorate, sodium 
chloride and sulphuric acid as raw materials, is out- 
lined. Recent improvements in the Mathieson, Solvay 
and Day-Kesting processes are reviewed, as well as 
the new Columbia-Southern modifications of the 
Persson process. 

Trends in absorption tower design, in materials 
of construction and in bleaching with chlorine dioxide 
are briefly indicated. 


HLORINE dioxide has come into wide-spread use as a 

bleaching agent in the pulp and paper industry in the 
last five years. It is explosive spontaneously at high par- 
tial pressures and readily at relatively low partial pressures 
by an electric spark, or any other means of producing a 
high temperature locally. Therefore chlorine dioxide 
cannot be manufactured at a central location and shipped 
to the consumer economically. For small uses such as 
flour bleaching where a high cost can be tolerated, 
chlorine dioxide is being frozen as a hydrate at 16% 
concentration and shipped and stored under refrigera- 
tion. ") For use as an antiseptic, chlorine dioxide is being 
dissolved in a solution of sodium perborate, with which 
it forms a stable complex ). A 4% solution can be safely 
shipped, from w hich the chlorine dioxide is released by 
acidification at the point of use. 

For large scale use chlorine dioxide is manufactured in 
the consumer’s plant from sodium chlorate, sulphuric 
acid and a reducing agent, which is usually either sul- 
phur dioxide or methanol. There are now 50 chlorine 
dioxide plants in pulp mills on this continent, 32 using 


1Manuscript received September 24, 1958 

Professor of Chemical Engineering, University of Toronto, Toronto, Ont. 
Based on a paper presented at the Joint A.I.Ch.E.-C.1.C, Chemical 
Engineering Conference, Montreal, Que., April 20-23, 1958. 
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sulphur dioxide, 17 methanol and one hydrochloric acid 
as the reducing agent. 

The six different processes now in use for generating 
chlorine dioxide in pulp mills were compared by the 
author in 1954), Three of these use sulphur dioxide 
directly as the reducing agent. In the Rapson-Wayman 
process “) concentrated sodium chlorate solution is con- 
tinuously passed down a packed tower counter-current to 
a stream of dilute sulphur dioxide in air or other inert 
gas, no sulphuric acid being required. The Holst pro- 
cess “), widely used in Sweden, involves treating a sul- 
phuric acid solution of sodium chlorate with dilute sul- 
phur dioxide gas batchwise in a tank. The Mathieson 
process (°) continuously generates chlorine dioxide from 
a sulphuric acid solution of sodium chlorate with sulphur 
dioxide in a vigorously agitated tank. In the Persson pro- 
cess “) the reducing agent is chromic sulphate, which is 
continuously oxidized by the acidified chlorate to chromic 
acid, which is reduced to chromic sulphate by sulphur 
dioxide in a separate cycle. The Solvay process (5) in- 
volves continuously reducing sodium chlorate with meth- 
anol in sulphuric acid solution. The Day-Kesting pro- 
cess “) involves oxidizing sodium chloride to chlorate in 
electrolytic cells, and reducing the chlorate to chlorine 
dioxide with hydrochloric acid in two interlocking con- 
tinuous cycles. 

Since 1954 several of these processes have been im- 
proved in important respects. Considerable progress has 
been made in understanding the mechanism of the re- 
actions by which chlorine dioxide is formed in these pro- 
cesses, and this has led to improved operation of the 
existing plants and the development of a new process. 


The mechanism of formation of chlorine 
dioxide from chlorates 


The author has previously suggested that chlorine 
dioxide is formed from chlorates by the same_ basic 
mechanism in all processes “°), This concept is supported 
by certain facts. First, in all processes the chlorine dioxide 
is formed from chlorates only in strong acid solutions. 
It is very significant that the Holst, Rapson- -Wayman, 
Solvay, Mathieson and Persson processes all operate in a 
solution which is between 8 and 10 normal in sulphuric 
acid. In the case of the Rapson-Wayman process, the 
solution ends up about 9 normal in sulphuric acid, even 
though no sulphuric acid is added in the process, and no 
attempt is made to control the acidity. 
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The second significant fact is that chlorides are found 
in the generator solution in all processes, usually in a con- 
centration between 0.02 and 0.08 molar. In the Holst 
process which is operated batchwise, the chloride con- 
centration is almost vanishingly low at the beginning when 
the chlorate concentration is high, but it increases gra- 
dually as the chlorate concentration drops, until it is 
relativ ely high when the chlorate is exhausted. 


Therefore it is postulated that the primary reaction in 
all chlorine dioxide manufacture from chlorates is the 
reaction between chloric acid and hydrochloric acid to 
form chlorine dioxide and chlorine. The following re- 
action mechanism would account for this satisfactorily 
by a series of bimolecular reactions; 


HCIO, + HCl HCIO,+HCIO..... (1) 
HCIO, + HCIO, > 2CI0,+H,O .... .~ (2) 
HOO 4+ NO46O,4HO......-@ 


Summary reaction 2HCIO, + 2HCI > 2Cl0, + 
DADE 6c bie ee wnc dwn OB 


This mechanism accounts for the stoichiometry of 
reaction (4) very well. Reaction (3) is well known. The 
hypothesis receives additional support from the fact that 
in our laboratory, conditions have been found for con- 
verting one mole of sodium chlorite into two moles of 
chlorine dioxide by reaction with an excess of sodium 
chlorate in sulphuric acid solution as in reaction (2) 
However, this mechanism cannot be said to have been 
proven. An intensive investigation is under way attempting 
to elucidate the mechanism, but the over-all reaction (4) 
has been clearly demonstrated to occur. 


The R2 process 


Arising out of this hypothetical reaction mechanism, 
a new process for making chlorine dioxide has been 
developed (), tested in pilot plants, and patented (*), in 
which sodium chlorate and sodium chloride in sulphuric 
acid are very nearly quantitatively converted into chlorine 
dioxide and chlorine according to reaction (4). In this 
continuous process, for convenience called the R2 process, 
streams of concentrated sodium chloride solution, con- 
centrated sodium chlorate solution and concentrated sul- 
phuric acid are metered into a vigorously agitated reaction 
vessel (Figure 1). Air is blown into the reactor through 
porous plates over which draught tubes are placed to 


—» C102 
Cla 


+ 
Effluent 


H2S04 
NaHSQ4 


Cl 
ce 


Air 


Figure 1—The R2 chlorine dioxide generator. 
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promote circulation of the liquid. Due to the vigorous 
turbulence, the concentrated solutions are immediately 
diluted by the generator solution, which continuously 
overflows from the vessel. The air continuously removes 
the chlorine dioxide and chlorine formed in the solution 
by reaction (4). A steady state is reached in which the 
composition of the generating solution remains constant. 

The optimum conditions within the reactor are 0.1-0.2 
molar sodium chlorate, 9-10 normal sulphuric acid, and a 
low chloride concentration. The latter seeks its own level, 
depending on the acidity, the concentration of the chlor- 
rate, and the temperature. Usually the chloride concen- 
tration falls within the range 0. 02 to 0.08 molar. The 
solution is very nearly saturated with respect to sodium 
acid sulphate. 

In the R2 process the chlorate which is decomposed is 
converted into about 99% of the theoretical chlorine 
dioxide. The ratio of chlorine dioxide to chlorine formed 
is very close to 2 moles to 1, indicating that the reaction 
proceeds nearly theoretically according to reaction (4). 

The chlorine dioxide is absorbed in water in a packed 
tower preferentially to chlorine due to tenfold greater 
solubility. Although about 10% chlorine remains in the 
chlorine dioxide, very little chlorine dioxide leaves the 
tower with the chlorine, which may be dissolved in 
sodium hydroxide solution to form sodium hypochlorite, 
or used in some other way. 


Very little heat is produced in the R2 generator, other 
than that evolved by dilution of the concentrated sul- 
phuric acid. Therefore the generator can be of very 
simple design, since no heat exchange is required. 

Unfortunately when sodium chloride is used along 
with sodium chlorate, twice as much sulphuric acid is 
required per unit weight of chlorine dioxide to purge 
the extra sodium ion from the system while maintaining 
the required acidity. Therefore if the sulphuric acid and 
the sodium sulphate i in the liquid effluent from the reactor 
are thrown away, the R2 process does not make chlorine 
dioxide as cheaply as other processes. If the effluent is 
used to provide sodium sulphate make-up for the kraft 
pulping process, it is cheaper than other processes in 
locations where sodium sulphate is expensive and sul- 
phuric acid is cheap, such as in southern United States. 
However, if sodium acid sulphate is crystallized out of 
the effuent and used for kraft make-up, the excess sul- 
phuric acid is recovered for re-use. The R2 process then 
makes chlorine dioxide cheaper than any previous process 
from sodium chlorate. 

The R2 process has been assigned to Hooker Electro- 
chemical Company in the United States and Electric 
Reduction Co. of Canada Limited in Canada. 


Application to other processes 

It is postulated that the reactions taking place in the 
R2 generator are those which produce the chlorine 
dioxide in all other generators. The only difference lies 
in the fact that in the other processes a reducing agent 
such as sulphur dioxide, methanol or chromic acid is used 
to ae the formation of chlorine, or at least to mini- 
mize it. This is lage ge by reducing either hy po- 
aaa acid or chlorine to chloride, thus maintaining 
the smmannaihed of chloride necessary for formation of 
chlorine dioxide by reaction (4). 


HCIO + H,SO, = HC] + H,SO, ..... @) 
HClO + CH,OH = HCl + HCHO +H,0 . (6) 
HCiO + HCHO = HCi+ HCOOH .... (7) 


If, due to the escape of some chlorine or chloride 
from the generator, the chloride concentration tends to 
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drop below that required for maintaining the desired 
chlorine dioxide production rate, or if excess reducin 
agent is fed into the generator, some chlorate will be 
reduced to chloride providing the necessary reactant. 
HClO, + H,SO, = HCIO, + H,SO, . . . . (8) 
HCIO, + H,SO. HCIO + H,SO, .... (9) 
HCIO + H,SO, = HCl +H,SO, ... . . (10) 


The importance of chloride in all chlorine dioxide 
manufacturing processes has been clearly demonstrated 
by adding silver sulphate to precipitate chlorides in 
laboratory experiments with acidified chlorate and various 
reducing agents. In no case was any significant amount 
of chlorine dioxide formed. Even with sodium chlorite 
in acid solution no chlorine dioxide was formed in the 
presence of silver sulphate, indicating that chloride is 
essential for the formation of chlorine dioxide. 

If the postulated mechanism is correct, then two 
sources of inefficiency in most chlorine dioxide manufac- 
turing processes are the reduction of some chlorate to 
chloride, and the loss of chlorine from the generator. 
Since both the chloride over-flowing in the effluent solution 
and the chlorine leaving with the chlorine dioxide arise 
from reduction of some chlorate to chloride, the author 
has suggested that the addition of sodium chloride to 
the chlorate in all processes may improve the efficiency. 
A patent has been applied for on this procedure. Olin 
Mathieson Chemical Corporation has recommended that 
5-10% sodium chloride be added to the chlorate solution 
fed into the Mathieson generator.*) 

This concept of the mechanism of the reactions pro- 
ducing chlorine dioxide has helped greatly in under- 
standing the operation of chlorine dioxide generators. For 
example in the Mathieson generator it is often found that 
the rate of chlorine dioxide production is not immediately 
increased on increasing the rate of flow of sulphur dioxide 
into the generator. This is now readily understood, since 
the rate of chlorine dioxide production cannot be increased 
except by increasing the chloride concentration, when 
the temperature, acidity and chlorate concentration are 
maintained constant. Therefore the increased rate of pro- 
duction of chlorine dioxide desired cannot be accomp- 
lished until sufficient chloride has been formed by reduc- 
tion of chlorate with SO, to provide the necessary 
chloride concentration for the new level of production 
rate. 

If the SO, flow is stopped without decreasing the rate 
of air flow, chlorine dioxide continues to be produced at 
a slowly diminishing rate as the chloride concentration 
diminishes in the generator, At the same time the chlorine 
content of the gas leaving the generator increases. This 
shows that reaction (4) is taking place. 

If a generator is allowed to stand overnight without 
feeding in any chlorate or sulphuric acid or sulphur 
dioxide, the chloride level is very low when production 
of chlorine dioxide is resumed. Therefore on start-up 
chlorine dioxide is not produced at the desired rate as soon 
as sulphur dioxide is fed into the generator, but production 
begins at a very low rate and gradually increases up to the 
desired rate. This phenomenon is accounted for by the 
need to build up the chloride concentration to the 
required level before the desired production rate can be 
established. This start-up period can be eliminated by 
adding to the generator sodium chloride in sufficient 
amount to immediately establish the chloride concentra- 
tion at the required level. 


Improvements in the Mathieson process 


A study of the Mathieson process in a pilot plant in 
our laboratory showed that under certain conditions some 
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sulphur dioxide gets through the generator solution un- 
reacted. On reaching the gas phase above the generator 
this sulphur dioxide first reacts with the chlorine present 
in the gas producing sulphuric acid and hydrochloric 
acid. (4) If the amount of SO, escaping is more than 
enough to reduce all the chlorine to hydrochloric acid, 
the sulphur dioxide will then react w ith some chlorine 
dioxide. The hydrochloric and sulphuric acids appear in 
the chlorine dioxide solution. 


The scrubber 


This difficulty was overcome by installing a scrubber 
on top of the generator consisting of a small tower packed 
with Raschig rings (Figure 2). The scrubber is designed 


To Absorption Towers 
~ ClO, + Air + (Cl,) 


Gos Scrubber 


To stripper 
and waste or 
=—* recovery 


Primary 
Generator 


Secondary 
Generator 


Figure 2—The Mathieson ClO. generator, as modified by 
Electric Reduction Co. of Canada Limited. 


in such a way that the chlorate solution passes down the 
packing on its way into the generator while the gas 
passing upward through the packing is scrubbed by the 
chlorate solution. In this way no acid is carried from the 
generator into the absorption tower. 


The scrubber then serves the purpose of returning to 
the generator the hydrochloric acid produced by reaction 
between chlorine and sulphur dioxide, along with the sul- 
phuric acid formed in the same reaction. This chloride 
helps maintain the chloride concentration in the generator, 
making it unnecessary to reduce chlorate to chloride for 
this purpose. Therefore the efficiency of the process is 
improved by one or two per cent. 

The same phenomenon of sulphur dioxide escaping 
from the generator liquid has been found in all commer- 
cial installations so far tested. In fact all Mathieson gener- 
ators, when proper analytical methods are used, have been 
shown to produce chlorine dioxide free, or nearly free 
of chlorine, but containing some hydrochloric acid. There- 
fore for the newer Canadian installations, Electric Reduc- 
tion Co. of Canada has designed a scrubber which has 
brought about the expected slight improvement in effi- 
ciency. (5) 


The secondary generator 


For economic operation of the Mathieson process the 
concentration of chlorate in the generator must be main- 
tained between 10 and 20 gm. per litre, depending on the 
production rate. If the chlorate concentration is dropped 
below this level, the efficiency is significantly lowered, 
because the chloride concentration must be increased to 
maintain the production rate. Since the effluent from the 
generator normally is wasted, 20 gm. per litre chlorate 
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solution flowing out of the generator represents a loss of 
about 5% of the ingoing chlorate. To recover part of this 
chlorate as chlorine dioxide, Electric Reduction Co. of 
Canada has designed a small secondary generator into 
which the effluent from the main generator flows by 
gravity. By means of a single porous plate and one draught 
tube a stream of sulphur dioxide is introduced into the 
secondary generator, producing chlorine dioxide at some- 
what reduced efficiency, but almost stripping the solution 
of sodium chlorate (Figure 2). 


The gas from the secondary generator goes to the top 
of the primary generator to be mixed with the gas 
produced there. If there is any excess sulphur dioxide or 
extra chlorine in the gas from the secondary generator 
they will react with either chlorine or sulphur dioxide in 
the primary generator gas, and the hydrochloric and 
sulphuric acids so formed will be scrubbed out of the gas 
in the scrubber and returned to the system. Thus the 
combination of the secondary generator and the scrubber 
produces chlorine dioxide with the maximum possible 
efficiency. (5) 


Improvements in the Solvay process 


The Solvay chlorine dioxide generator has been 
completely re-designed since ‘the last comparison of 
processes was made. The three packed towers in which 
the sulphuric acid solution containing chlorate and 
methanol was circulated provided many operating prob- 
lems, especially the maintenance of the pumps circulating 
the extremely corrosive liquid.) They have now been 
replaced by two jacketted lead- lined steel reaction vessels 
operated in series in cascade fashion, over-flowing by 
gravity from one to the other (Figure 3). Each vessel is 


To Absorption Towers 
ClO, + Air +(Cl,) 


To waste or 
recovery 
Hot 


Figure 3—The new Solvay ClO, generator. 


agitated vigorously by an air stream, which also serves 
to remove the chlorine dioxide from the solution and 
carry it to the absorption tower at low partial pressure 
to prevent explosion. 

Concentrated aqueous sodium chlorate solution, 
methanol, and concentrated sulphuric acid are metered 
into the first reactor and immediately mixed with the 
generator liquor due to the vigorous agitation. Chlorine 
dioxide is produced at a constant rate in this solution 
which is maintained in a steady state at about 100 g.p.l. 
sodium chlorate and 9 normal sulphuric acid. The solution 
overflows to the second generator which is maintained 
in a different steady state, with lower chlorate concentra- 
tion, perhaps 10 to 20 g.p.l., by adding methanol con- 
tinuously. 
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Since most of the reaction takes place in the first 
generator, it is necessary to circulate cool water through 
the jacket to maintain about 140°F. Since the smaller 
portion of the reaction takes place in the second generator, 
it is usually necessary to add heat to the second generator 
to keep it about 145°F. 

This new design of the Solvay process has improved 
the efficiency of its oper ation. With careful control of all 

variables at their optimum values it is possible to obtain 
90% efficiency. 


The Day-Kesting process 


A Day-Kesting plant integrating manufacture of 
sodium chlorate and chlorine dioxide has been built at 
Berlin, N.H., for Brown Company, incorporating several 
improvements over the similar European plants.*) The 
chief improvement is the use of one large common 
storage tank from which the sodium chlorate solution is 
circulated rapidly through the electrolytic cells for oxi- 
dation of chloride to chlorate (Figure, 4). From the same 
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2NeC/ + 640 —_ 2NOCI0y + Sy 2 NOCI0y + 4HCl—-~ 2NOCI + 2C1Op+C lg s2ryO 
Figure 4—The Day-Kesting ClO. process. 


tank the solution is also circulated at a very slow rate 
through the six ClO, reactors arranged in cascade fashion. 
A large part of the chlorate in the solution is reduced to 
chloride by reaction with HCl forming chlorine dioxide 
and chlorine, and the partially exhausted sodium chlorate 
solution is returned to the storage tank and mixed with 
the main solution. 


The electrolytic cells are about two-thirds larger and 
the chlorine dioxide reactors three times larger than those 
in Europe, decreasing the number required to 20 cells and 
six reactors for the production of 13 tons of ClO, per day. 


This process has not been widely adopted because of 
the high capital investment required. Although the chlo- 
rine dioxide is cheaper than that made from purchased 
chlorate, the pay-out time required to recover the original 

capital is longer than pulp companies are apparently 
willing to accept. 


The Columbia-Southern process 

The Persson process has been investigated thoroughly 
by Columbia Southern Chemical Corporation and re- 
designed to make a new process with some novel features 
(Figure 5). It is proposed to remove the chlorine 
dioxide from the generating liquor by steam instead of 
air. At first sight this does not seem feasible, since chlorine 
dioxide normally decomposes spontaneously at 100°C. 
However, the concentration of chlorine dioxide in the 
steam is kept low enough to prevent this decomposition. 
It is therefore suggested that the chlorine dioxide can be 
used directly with the pulp instead of dissolving the 
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Figure 5—The Columbia-Southern ClO: process. 


chlorine dioxide in water solution before mixing with 
the pulp. If the steam and chlorine dioxide can be added 
to the pulp simultaneously in this manner it would avoid 
the dilution of the stock by a large amount of cold water. 
However, abandoning the storage tank would involve 
some operating problems, since the chlorine dioxide pro- 
duction would have to be very accurately balanced with 
the pulp production, and there would be no surge capacity 
to make up for any lost time in the chlorine dioxide plant. 

The Columbia Southern process uses a single generator 
instead of the four normally provided in the Persson 


process. An excess of trivalent chromium is used to - 


convert all the chlorate in a single pass through the 
generator, but this creates some difficulties in the crystal- 
lizing system, which is an essential part of the process. 


Absorption towers 


Many chlorine dioxide plants have been built with 
absorption towers of too large diameter to operate effi- 
ciently at the water and gas flow rates required for the 
normal production rate of the generator. Since packed 
absorption towers have a relatively narrow optimum 
operating range in this application, significant amounts of 
chlorine dioxide escape from some towers. This loss is 
increased sharply when the plant is operated at less than 
its maximum rate. One 24 in. diameter tower is adequate 
for four tons ClO, per day but at one ton per day it is 
too big. 

Therefore two absorption towers of different diameter 
should be provided to give three optimum operating 
rates, with one or the other tower, or with both in parallel. 
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Materials of construction 

The experience with respect to materials of construc- 
tion for chlorine dioxide production and bleaching in 27 
pulp mills has been reported recently by W. H. Rapson 
and R. J. Neale “’). Homogeneously lead-lined steel has 
been most widely used for chlorine dioxide generators, 
but a recent installation of glass-lined steel appears to be 
satisfactory. For piping, glass, glass-lined steel, saran-lined 
steel and unplasticized polyvinyl chloride have been found 
satisfactory. For pumps Durichlor, Worthite and Ircamet 
are w idely used. Saunders-type valves with saran-lined 
or glass- -lined steel bodies and Teflon or Kel-F diaphragms 


are most common. 


The biggest corrosion problem in using chlorine di- 
oxide has been the equipment for mixing it with the hot 
pulp. Although titanium is expensive, it is the only 
material that seems to give long life in this application. 
A recent development has been: stainless steel mixers 
coated with fibreglass and polyester resin. 


Bleaching with chlorine dioxide 


Improvements in equipment and processes which make 
chlorine dioxide more cheaply are very important to the 
pulp and paper industry. This bleaching agent is unique 
in its ability to bleach pulp very white without damaging 
the strength of paper made from the pulp, and without 
introducing carbonyl groups into the cellulose molecules 
which cause the white pulp to become yellow with 
age (19, 20), 

The modern trend in bleaching is to eliminate hypo- 
chlorite altogether and bleach in five stages with chlorine, 
caustic soda, chlorine dioxide, caustic soda and chlorine 
dioxide. These five treatments, applied under suitable con- 
ditions for each, bleach kraft pulps, soda pulps, sulphite 
pulps and neutral sulphite pulps to a new level of white- 
ness, color stability, strength and uniformity which has 
created a higher standard for pulps in the paper industry. 
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Catalytic Reforming’ 


N. J. EMMS? 


This paper describes the Catalytic Reforming 
Process for upgrading naphtha feedstocks to high 
octanes using a regenerable platinum catalyst. Par- 
ticular reference is made to the Powerforming 
Process developed by Esso Research and Engineering 
Company and in use by Imperial Oil Limited. 

The various types of fixed bed units are described 
and a simplified flow plan of a typical Powerforming 
Unit is illustrated. The general operation of the unit 
is discussed as well as the primary reactions involved 
in catalytic reforming. Operating conditions are listed 
and some product yields and quality are shown in a 
table for various feedstocks. 

Regeneration is required to reactivate the catalyst 
by burning off carbon and sulphur compounds laid 
down during the oil cycle. A general description of 
how this is carried out on a unit containing a swing 
reactor is outlined. 


i the past few years octane requirements for auto- 
mobile engines have risen rapidly due to the use of 
higher and higher compression ratios by the automobile 
manufacturers. In order to meet this demand for high 
octane gasolines it has been necessary for the petroleum 


TABLE 


industry to develop improved refining processes. Processes 
like catalytic cracking, which were capable of making 
satisfactory high octane gasoline blending stocks just a 
short time ago, are incapable of making the higher octanes 
required. One of the most recent developments to come 
into widespread use is the Catalytic Reforming Process. 


The Catalytic Reforming Process is a continuous 
catalytic process to upgrade low octane virgin or thermal 
naphthas, which can no longer be tolerated in gasoline 
blends, into high octane components. As a result of the 
intense competition in the oil industry, there were as 
many as 15 or more individual types of processes in use 
as of last year. The majority of ‘these designs are of the 
fixed bed type but the use of the fluid and moving bed 
principles have been used in a few designs. It w ould be 
impossible to review the individual designs in a paper 
such as this and therefore the various processes are listed 
for record purpose in Table 1. 


At Imperial Oil Limited, the Powerforming process 
developed by Esso Research and Engineering Company 
has been adopted and units installed at nearly all their 
refineries. In the following discussion, therefore, on 
catalytic reforming, the Pow erforming process will be the 
basis of the comments. 


1 


CATALYTIC REFORMING PROCESSES 


Process Developed By Type 

Cat ORIG. oc ene ciines 38 DeRo eased ee ed PRAMANNE DE CERISE OAR ooo iaho ere Sad x sicls ered deme d ea aoe .Fixed Bed 
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PRON ic ovine coasts naEweosure ees RFI RONUGE CE PINMIIE UB iiss oie odes cso nG cee Ceca ten heeuaeas . Fixed Bed 
NTI SE PIO 56 od ok Siale Ss » vincw'w 00 wei o oe < RINE EEO she aft iracinl oe dhes- see MERE hard a ae asennad ee a . Fixed Bed 
SOVAIOEMHNG «566i 5 io sors hae totais ce ees + vCORY ROOM EEE AMES MEM ESE RS 2 Fixed Bed 
Fixed-Bed Hydroforming................. Esso Research & E ngineering Company and M. W. Kellogg Co...... Fixed Bed 


Pow erforming. . 
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Pirate daca etre tate Esso Research & Engineering Company.....................+06. 
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.Fixed Bed 
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Fluid sydeekenming Aa aial n wesrete chattel ieee tea Esso Research & Engineering Company..... i cwesee 4 cieoe a nae 
INI sivas own edad w alee nee Oe ee 1 AC CT eee eee ene ea rer eee err eee Fluid Bed 
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Based on a paper presented at the 41st Annual Conference and Exhi- 
bition, The Chemical Institute of Canada, Toronto, Ont., May 1958. 
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Types of fixed bed reforming units 

In a fixed bed catalytic reforming unit, there can be 
three classifications. The first could be called a non- 
regenerative unit which is designed to operate for a year 
or :nore on the original catalyst charge. When the catalyst 
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becomes inactive it is removed from the reactors and sent 
to the catalyst manufacturer for restoration to its original 
activity level. The second type could be called a semi- 
regenerative type. This is essentially the same as the non- 
regenerative type but has added facilities for the regenera- 
tion of the catalyst at about a three month interval. 
Regeneration of the catalyst consists of burning off the 
carbon which has been laid down on the catalyst and 
resulted in the deactivation. In the semi-regenerative unit, 
feed must be removed from the unit during regeneration. 
The third type of fixed bed, called the regenerative type, 
is based on frequent regeneration and continuous operation. 
A swing reactor is included in the design to allow indi- 
vidual reactors to be taken off stream for re generation at 
which time its place in the sequence is taken by the 
swing reactor. 


Process flow 
A typical fixed bed regenerative Powerforming unit is 


shown in Figure 1. This is a much simplified flow 
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Figure 1—Simplified powerforming unit flow plan. 


diagram for easy reference. The process flow through this 
unit is as follows. The first step in the process is the 
preparation of the naphtha feed. Fractionation equipment 
often precedes the unit and in this way the correct boiling 
fraction is distilled for charging to the unit. Since sulphur 
is a poison to platinum type c: ataly sts, the naphtha fraction 
is normally treated to remove sulphur compounds. As 
shown in the process flow diagram desulphurization is 
carried out in a hydrofiner w hich consists of a fixed bed 
reactor filled with a cobalt molybdate catalyst. The 
naphtha feed, heatcd to around 550-650°F., is mixed with 
fresh hydrogen-containing gas from the reforming unit 
and pumped through the reactor. After leaving the 
reactor the product. is cooled and taken to a gas- liquid 
s:parator for product and gas separation. The hydrogen 
rich off-gas can be used other hydrogen processes or 
put to fuel gas. 

After separation the liquid product is pumped through 
an absorber stripper tower where any residual hydrogen 
sulphide is stripped from the liquid. The liquid product 
then combines with recycle gas and passes through the 
preheat furnace and the first reactor. In this flow plan four 
reactors are shown as being onstream with a spare reactor 
as a swing for regeneration, Between each reactor is a 
furnace for rehe: ating the feed to the required temperature 
before charging to each reactor. From the last reactor the 
product is cooled and sent to a product separator for 
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liquid and gas separation. The liquid product from the 
separator is then charged to a debutanizer to stabilize it 
before it is sent to storage or blending. Part of the 
hydrogen-rich gas from the separator is compressed and 
recycled to join with the feed to the preheat furnace and 
first reactor. 


Primary catalytic reforming reactions 

At this time it is probably opportune to comment 
briefly on the primary reactions involved in catalytic 
reforming. The basic idea, of course, in reforming is 
the conversion of low octane naphthenes and paraffins to 
high octane aromatics. The higher the aromatic content 
of the product, the higher the octane. In obtaining this 
conversion a yield loss is involved. The first three reactions 
shown in Figure 2 are the primary processes for octane 
improvement. 


(1) Naphthene Dehydrogenat lon 


CH- 
2 
Cyclohexane Benzene 


This reaction is endothermic and rapid 


2 Naphthene Dehydroisomerization CH, 
4 
CH 


3 Cc. 
C“H ” CH; . 


Ng 
CH 

Ha 

1, 2 dimethylcyclopentane methylcyclohexane toluene 


This overall reaction is endothermic 


3 


3) Paraffin ae” Cy 


sinc lll mene 
CH 3(CHip) 5 CH3 =F 
CH 


n-heptane methylcy clohexane toluene 


This reaction is strongly endothermic. 


Figure 2—Primary catalytic reforming reactions, 


Other reactions that take place are the isomerization 
of paraffins, hydrocracking and the saturation of olefins. 
Although these latter reactions are exothermic the main 
reactions are endothermic which accounts for the drop 
in temperature occurring across the reforming reactors 
and the need for the reheat furnaces. 


Typical operating conditions and catalysts 

Typical operating conditions on a Powerforming unit 
are 300-600 psig pressure, 900-1000°F. reactor inlet tem- 
perature, a 4000-6000 SCF/bbl. recycle gas rate and 1-4 
space velocity. (Pounds per hour of product per pound 
catalyst ) The heart of the process can probably be called 
the cataly st. A number of catalysts have been developed 
by the various licensors, with composition of the catalysts 
varying but the majority in use today are platinum- 
containing catalysts. One of the first reforming catalysts 
was a mixture of molybdenum oxide and alumina, and a 
fluid molybdena and alumina catalyst is used in the fluid 
hydroforming units. In the moving bed Hyperforming 
process a cobalt molybdate pellet on a silica stabilized 
alumina base is used while in the Thermofor Catalytic 
Reforming Process a synthetic bead of co- precipitated 
chromia and alumina is the c catalyst. The majority of the 
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fixed bed, continuous catalytic reforming processes use 
an acid type, pelleted cataly st containing platinum (0.01 
to 1.0%) aver aging approximately 0.5% on alumina or 
silica alumina carriers. 

The catalyst used in Powerforming was developed by 
the Isso Research and E ngineering Company and consists 
of platinum on an alumina base. The catalyst has a high 
activity and a long life. A platinum c catalyst can be 
permanently poisoned by heavy metals such as molyb- 
denum, arsenic and lead. This is usually not a problem, 
as the feeds normally prepared are overhead distillates 
from a fractionation tower. Sulphur is a temporary 
platinum catalyst poison as well as an impurity that can 
lead to a severe corrosion of process equipment. It is for 
this reason facilities such as the hydrofiner in the process 
flow diagram are normally installed for desulphurization 
of the feed. Any scale formation from corrosion tends to 
build up at the top of the reactors leading to pressure 
drop problems. 


Regeneration 

At this point a few comments on regeneration would 
be worthwhile as they apply to a fixed bed unit of the 
regenerative type where a swing reactor has been included 
in the design to allow for a continuous operation. Catalyst 
regeneration is the name applied to the technique of 
burning off the carbon and sulphur compounds that have 
been laid down during the oil cycle and deactivated the 
catalyst. Regeneration essentially restores original catalyst 
activity and selectivity. The reactor to be regenerated is 
replaced in the reaction sequence by the swing reactor. In 
such a process as reforming where a large volume of 





hydrogen gas is being circulated special safety controls 
are required. To prevent the hazard of mixing regenera- 
tion gases and hydrocarbon vapors, the reactors are often 
provided with double block valves and vents. This is 
provided on the Powerformers of this type built for 
Imperial Oil Limited. The double valve system is of great 
importance as it makes possible the blanking off of one 
reactor from service for repairs or ¢: atalyst change without 
affecting the rest of the unit. The hy drocarbon valves for 
each reactor are normally motor operated and push 
button controlled from the control One of the 
regeneration values is motor operated and the other is 
hand operated. These valves are interlocked electrically 
to prevent errors in the sequence of valve operation. 
Oxygen recorders are required on the inlet and outlet 
regeneration gas stream and each recorded is provided 
with temperature recorders for observing regeneration 
tsmperatures. Sometimes a seal gas system is provided 
which uses inert gas to maintain an elevated pressure 
between the double block valves to insure isolation of the 
reactor to be regenerated from the on-stream equipment. 


room. 


The reactor to be regenerated, after removal from the 
sequence, is purged of heavy hydrocarbons by depressur- 
ing and pressuring with inert gas. After purging, a flue 
gas containing approximately 2% of oxygen is passed 
over the cataly st at a temperature between 800°F. and 
1050°F. After the coke burn is completed, the reactor is 
purged with inert flue gas followed by a recycle gas 
purge and is placed back in reaction. The period required 
for regeneration from the time a reactor is taken off 
stream to the time it is put back will vary but experience 
has shown this sequence can be carried out in 16-24 hours 
if necessary. 


TABLE 2 


TYPICAL 


Feed Source Arabian 
INCH E UG RRO lo ssiers. Haid < Olmmdo o Rae Be 160/260 
Feed Properties 
Gravity °API. 64.7 
Research Oc tane Cle ir Ce a ee eee ‘ 44 
AStM S% Point “F....:.. Sb hewa a aaNet ania ks 5 199 
PAU te ee MORNE Beats Sy cid s Sie btw boehladicieeie ln wire ns : 217 
ASTM 95% Point °F..... Scr iat Dieta Sh ait 249 
RCRD Mee OHI) Oy aia sy ct a ss pe eee wee 32:2 
Yields ©% on Feed 
Meee We ts oh Seto a aes be gta eh Rte ee ot ee 60.1 
BN WE aria tova'y suet oe) Liv Siho Ne Sra eS emo er eet 1.6 
MMR cc A ae eS ie ie ee eS 2.9 
C2, Wt. % 7.0 
C3 Wt. % 11.0 
C4, Vol. % 17.0 
C5+ Powerformate Properties 
Volatility, © Distillation plus Loss @ 158°F. IK 
ais: $3 
257-6... 77 
6 al ere 97 
356°F. 98 
RMVIEW He 8 wh Sek ches. Capa deu imeem. 50.7 
Reid Vapour Pressure p.s.i. 6.8 
Octane Numbers 
Research Clear... er oy 100 
+ 3cc Tetra ethyl lead. eteics 105.8 
Motor Clear. 88.6 
+ 3cc Tetra ethyl lead. 94.8 
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YIELDS AND QuaLity Data 


Mixed Louisiana Leduc 


Redwater 


170/310 200/330 200 ‘400 
52.2 55.0 54.8 
58 58 50 

220 223 247 

241 253 276 

288 304 341 
11.9 11.8 11.9 
tao 1 77.4 
2.4 nal 2.4 
2.1 1.7 1.8 
Sat 2.9 3.0 
6.2 5.3 4.9 
10.7 8.8 $7 
8 9 7 
27 29 7} 
Sa:5 50.5 36 
88.5 83 68 
98 Q7 O4 
$5.5 14.6 5.4 

3.9 }.? 3.0 

100 100 100 
105.0 105.0 105.0 
88.6 88.6 88.6 
93.6 93.1 93.6 
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Yield and product inspections 

Typical feedstocks for catalytic reforming consist of 
low octane naphthas usually in the boiling range of 175°F. 
to 400°F. As mentioned previously prefractionation equip- 
ment is often installed to make the desired distillation cut 
for the reforming unit. The yield of reformed gasoline 

varies depending on the severity of the reforming opera- 

tion. The severity of operation will depend on the par- 
ticular refinery situation but current Powerformers are 
being designed to produce a reformate of 100 Research 
Octane Number clear. Table 2 indicates some typical 
product yields and quality. 

At times it may be desirable to split a wide-boiling- 
range naphtha into separate fractions which can be 
reformed under different conditions. A high final boiling 
point results when reforming a high end point feed. This 
product would be rerun before inclusion in the gasoline 
pool. 


Future processes 

In conclusion, reference is made to the two processes 
listed in Table 1 as combination processes. These designs 
use the reforming unit as the start of the design and add 
to it another step. The first of these is known as as Rex- 
forming and was developed by Universal Oil Products. 
In this design an aromatic extraction process has been 


added to a conventional Platformer. The reforming section 
of the unit is operated at lower severity and the reformate 
is charged to an extraction unit for the removal of aro- 
matics. The low octane raffinate is recycled back to the 
reforming unit. A glycol solvent is used in the extraction 
section. This process can be modified further to bring out 
aromatics for chemical use. 


The second combination process is known as Iso-Plus 
Houdriforming and is licensed by the Houdry Process 
Corporation. The first part of this process is a conven- 
tional Houdriformer operated at moderate severity in 
conjunction with one of three possible alternates. The first 
of these is similar to the Rexforming Process where the 
reformate is extracted and the raffinate recycled back to 
the reformer. In the second case the reformate is extracted 
and the raffinate reformed separately. The third alternate 
takes the reformate and charges it to a thermal reformer 
with the C3 and C4 olefins from thermal reforming being 
charged to a catalytic polymerization plant. 


The combination process is a direct result of the need 
for producing higher and higher octane components for 
motor gasoline. It is probable that for octane numbers 
substantially above 100 Research Octane Number Clear 
more processes along this line will be developed. 


~~ * 
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The Manufacture of Chemical 
Cellulose from Wood’ 


M. WAYMAN? 


This paper discusses progress made and problems 
associated with the manufacture of chemically pure 
cellulose from wood. The yield of pure cellulose 
which can be obtained from any wood species is 
limited by the chemical composition of the wood 
starting material. An examination of the chemical 
composition of various native wood species emphasizes 
the special potential of the genus Populus. Improved 
yield and good economy of manufacture of cellulose 
depends on the development of purification processes 
which are effective on less thoroughly digested wood. 
Conditions are developed for the most effective way 
of carrying out extraction of raw pulp with hot caustic 
soda. These include very high consistencies and very 
short times at high temperatures. The best conditions 
for extraction with strong caustic soda at low tem- 
perature are also reviewed. The need for properly 
designed equipment for both these processes is 
stressed. The special problem of removal of mannan 
is discussed. Two objectives are put forth: a more 
economical cellulose with today’s standard of purity 
and chemically pure cellulose from wood. 


ELLULOSE is a chemical, a long chain anhydro polymer 
oo glucose, which occurs in nature in fibrous form. 
The long fibres of cotton are nearly pure cellulose, and 
are readily spun into yarn, and the yarn woven into 
fabrics (1. 2), The linters associated with the cotton seed 
must be purified for use for chemical purposes. There 
are several processes available ‘%), 

Cellulose is also the main structural material in wood, 
and most cellulose for chemical uses is derived from 
wood (+ 5), The processes of purification of cellulose 
from wood form the substance of this paper. 


Purified cellulose, whether from wood or cotton, is an 
important raw material for several industries. It may be 
dissolved, usually by conversion to a derivative, and the 
cellulose regener rated from solution in many useful shapes, 
as fibres (rayon), film (Cellophane), sponges, ribbons, 
etc. Many of the derivatives are themselves useful. Cel- 
lulose acetate fibres are the most beautiful of the artificial 
silks. Fortisan, a cellulose regenerated from the acetate, 
is the strongest fibre. Cellulose acetate is also made into 
film, sheet, extruded and molded plastics. Molding pow- 
ders are often made of other esters, such as the propionate, 
or mixed esters, such as the acetate-butyrate. Cellulose 
ethers, methyl and ethyl cellulose, have found usefulness, 
as has carboxymethyl- cellulose, which is used in soaps to 


1Manuscript received August 24, 1958. 
2Technical Director, Columbia Cellulose Co. Limited, Prince Rupert, B.C. 
Based on a paper presented at the Joint A.I.Ch.E.-C.I.C. Chemical Engin- 
eering Conference, Montreal, Que., April 20-23, 1958. 
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retard soil deposition and in ice cream and in oil-well 
drilling mud as a thickener. Cellulose fibre properly puri- 
fied is often incorporated into urea-formaldehyde and 
melamine-formaldehyde moldings, and into artificial 
leathers. : 

Cellulose nitrate is used as a film-former in lacquers, 
as an ingredient of smokeless powder, and as a propellent 
in rockets. 


Wood as a raw material for cellulose 


This paper discusses progress made and problems as- 
sociated with the manufacture of chemically pure cellulose 
from wood. It will emphasize those steps in manufacture 
concerned with raising the purity to high levels, treating 
the pulping of wood and the bleaching of pulp only to the 
extent required to make a coh:rent picture. 

Cellulose is a carbohydrate. The total carbohy ae 
content of wood is about 70 to 80%, but of this the 
tion truly useful as a chemical is called the “alpha an 
lose”, which is generally only about 45% of the wood. 
Thus to purify cellulose from wood, more than 50% of 
the wood substance must be dissolved aw ay. 

The purification is generally accomplished in several 
steps. The wood must first be reduced to individual 
fibres. This is done by a process of digestion. There 
are two processes in common use for cellulose manu- 
facture, the sulphite and sulphate. By these processes, 
most of the undesired components of the wood are dis- 
solved, allowing the individual fibres to be separated. The 
individual fibres still contain unwanted impurities, mainly 
lignin and short-chain or odd-sugar carbohydrates. They 
are treated with chlorine to remove the lignin, and with 
caustic soda to remove the non-cellulosic carbohydrates. 
The cellulose so obtained is then bleached to make a white 
product. 

The amount and quality of cellulose obtainable from 
wood varies with the different species. Timell has re- 
cently given the cellulose content of various wood 
species (), and from his work and our own (‘), Table 1 
has been constructed. 


Several conclusions may be drawn from Table 1. 

1. It is at once evident that, except for the genus 
Populus, one cannot hope to get more than 44.5 to 
49.6% of alpha cellulose from wood; and that if 
rigorous exclusion of nonglucan material is ne- 
cessary, as it most certainly is in some applications, 
then these figures fall to 41-47%. These figures 
assume a complete separation of ‘es desired from 
the undesired constituents of wood in perfect 
yields. 

2. It is also evident that the genus Populus occupies 

a special place, in that its cellulose content is 2.4 

to 12% higher than the other woods (calculated 
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TABLE 1 


\LPHA CELLULOSE From Woop 


Alpha 
Species \lpha Cellulose 
Cellulose Corrected: 


‘ 
( 


Mannan) Xylan 


Aspen* 


Populus tremuloides 


Cottonwood 
Populus trichoca) pa 


White birch* 
Betula papyrifera 


Balsam fir 
Abie Ss balsame a” 
Abies amabilis 


Eastern hemlock* 
Tsuga canadensis 


Western hemlock 
Tsuga hete roph ylla 


Jack pine* 
Pinus banksiana 


White spruce* 
Picea glauca 


Sitka Spruce 
Picea sitchensis 


*From Timell 
: Corrected for non-glucan material 


on wood or 5 to 20% higher than softwood cel- 


lulose). The use of aspen, poplar and cottonwood 


is becoming quite widespread in the cellulose in- 


dustry. 

Use of cottonwood in the cellulose industry carries 
with it a number of problems — the high xylan 
content and short fibres characteristic of deciduous 
woods, and a resin with unusual chemistry, but it 
seems quite certain from both the laboratory and 
the plant that the extra cellulose yield is there and 


can be obtained in useful form. 


The table also shows that the deciduous woods are 
quite high in xylan content compared with the 
conifers. This has led to some special methods 
of treatment. 


On the other hand the mannan content of the 
deciduous woods is much less than that of the soft- 
woods. Mannan content can be quite critical 
especially in cellulose which is to be converted 
into cellulose acetate’®: !%, and sometimes extra- 
ordinary measures have to be taken to get rid of it. 


The deciduous woods generally do not pose this 
problem. 


Digestion 

In addition to cellulose, wood contains a high per- 
— of undesired components. These are lignin (16 

» 24, in hardwoods, 27 to 30% in softwoods), the non- 
an carbohydrate constituents (17-19%) and solvent- 
soluble extractives. These must be removed in making 
cellulose from wood. To remove the unwanted com- 
ponents, the physical form of the wood must be altered 
to convert it to a pulp of individual fibres. This is done 
by dissolving out nearly all of the lignin which cements 
the fibres together, in a process of digestion with acid or 
alkaline liquors, as explained earlier. 
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In this first step in cellulose manufacture, the bark 
is removed mechanically from the wood, which is then 
cut up into chips of such a size as to retain the fibre 
structure, and yet to make all the wood readily accessible 


to chemical treatment. Such chips are usually 4 in. thick 


and about & in. long, the other dimension being not 
critical. ‘The chips are then fed into pressure vessels, now 
commonly up to 19 ft. in diameter and 60 ft. high, and 
often stainless steel lined. Delignification is commonly 
accomplished by treating with an acid solution of a bi- 
sulphite, usually calcium, magnesium, sodium or am- 
monium(!), at about 150°C. 

Heating in the cooking operation is carefully con- 
trolled lest the reaction proceed too qnickly to com- 
pletion on the outside of the chips before the liquor has 
had a chance to penetrate to the inside of the chips. In 
this country four hours is considered a minimum in going 
from 70°C. to 150°C. A few plants elsewhere are oper- 
ating below that minimum. 

Delignification is largely complete shortly after reach- 
ing maximum temperature. The reaction appears to pro- 
ceed in two stages, a sulphonation of the lignin followed 
by a hydrolysis of the carbohydrate material. The first 
reaction is essentially complete when penetration is com- 
plete at 110 to 120 ‘C., and the higher temperature is re- 
quired for the carbohydrate hy droly sis. 


These reactions may be very fast, and the problem of 
their accomplishment is in part a physical one of pene- 
trating the wood structure with properly buffered treat- 
ment liquor and with heat, in that order. The liquor is 
not constant in composition, having an excess of SO: 
in solution, and having basic components which are taken 
up by the lignin sulphonic acids and wood acids as they 
are formed. The wood contains “air bubbles” which are 
hard to get rid In fact, normal coniferous wood at 
50% moisture content is more than half air. Many in- 
genious schemes have been developed to solve these pene- 
tration problems !?), 

Anyone familiar with the important role of wetting 
agents in the progress of the textile industry over the 
past 25 years or so would be tempted to suggested pene- 
trants. Some have indeed been tried, and a very modest 
degree of improvement has been achieved with one or 
two(!3), 


Wood may also be delignified by digestion with alkali, 
or more usually in an alkaline solution of sodium sulphide. 
This process, called the “sulphate” process, because so- 
dium sulphate is the purchased raw material for making 
the sulphide, has economic advantages in that the lignin 
is burned in a recovery system which is a source of cook- 
ing chemicals and heat. However, for making cellulose 
of sufficient purity for chemical conversion, the process 
is often modified by preceeding it with a “prehydrolysis” 
stage. This is heating in a solution of weak acid, either 
added mineral acids or simply organic acids formed from 
the wood, carried out in the same digester. The prehy- 
drolysis attacks the xylan to a far greater degree than the 
subsequent alkaline stage, making this process particularly 
useful with hardwoods. 


There are other chemical methods of delignification, 
such as digestion with acetic acid and several of its de- 
rivatives, and with alcohols to which some acid has been 
added. Some processes take place at atmospheric pres- 
sure, such as those with chlorite, chlorate or chlorine 
dioxide, peracetic acid and strong nitric acid. None of 
these has yet been developed to the point of commercial 
significance. 

In addition to the chemical and mechanical consider- 
ations associated with delignification, there are economic 
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considerations. Processes which are rapid and continuous 
shou!d require smaller capital investment than the present 
slow, bulky batch processes. Possibly the most important 
economic factor is yield. The percentage of pure cel- 
lulose which can be obtained from a tree is limited by 
the chemical composition of the wood as shown above. 
If the delignification process is carried on so far that 
cellulose hydrolysis occurs, and the amount of fibre re- 
covered from the digester is already less than the original 
cellulose percentage, and that in turn is perhaps only 85% 

pure, then nothing that is done subsequently can make 
up for what has already been lost. It therefore is im- 
portant to develop purification processes which are ef- 
fective on less thoroughly digested wood. This is the 
only way that theoretical yields and good economy can 
be approached. 


Chlorination 

After delignification, raw wood cellulose has a chem- 
ical composition which does not yet equal raw cotton 
linters. It has a lignin content of 1.5 to 3%, depending 
on the degree of delignification, 3 to 7% of mannan plus 
xy one and 0.5 to 1% of solvent-soluble extractives, as well 

2 to 4% of glucose-derived material which is not alpha 
ilee The alpha cellulose content is generally 86-98% 
(Table 2). 

TABLE 2 


TypicaAL COMPOSITION OF RAW CELLULOSE 


Raw Cotton?) | Raw Wood Pulp 
€ 


( ( 

Cellulose 94.0 i 86-89, 

Lignin _ 1.5-3 

Protein 1.3 

Mannan 4.0, 
2.9, 

Xylan 1.8, 
4.1, 

Pectic substances 0.9 

Solvent-solubles extractives 1. 4a 0.5-1.0 


refers to reference number 2. 
a alpha cellulose. 
» softwoods (of B.C. North Coast). 
« hardwood (cottonwood ). 
4 Wax plus organic acids. 


After blending and screening to remove improperly 
digested chips, knots and undissociated fibre bundles, the 
purification of cellulose from raw wood pulp generally 
follows a pattern: chlorination, caustic soda extraction and 
bleaching. The action of chlorine is to make the residual 
lignin water- or alkali-soluble. The caustic soda extrac- 
tion has the function of removing the short-chain or odd- 
sugar non-cellulosic carbohydrates, as well as the last 
traces of chlorinated lignin. The bleaching is, of course, 
to whiten the product, and indirectly the conversion 
products, that is, yarn, plastics or whatever is made from 
the cellulose. Some variants on this theme will be ex- 
plored. 


In order to remove the rest of the lignin it is usual 
to treat the woodpulp in dilute aqueous suspension with 
gaseous chlorine, or*with an aqueous solution of chlorine. 
The woodpulp and chlorine are mixed and pumped to the 
bottom of an up-flow tower, which overflows to a washer. 
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The operation is carried on continuously and at the 
prevailing water temperature, and should be short and 
at low pH. Any tendency to prolong the treatment, or 
any rise in pH at this point will cause oxidation rather 
than chlorination. Such oxidation affects adversely the 
properties of cellulose acetate solutions, and may make 
the final cellulose unsuitable as a starting material for 
cellulose derivatives. 


Caustic soda extraction 


The chlorinated pulp is well washed to get rid of a 
large part of the chlorinated lignin, and is then extracted 
with caustic soda solution. This removes the rest of the 
lignin, and, if properly carried out, nearly all the other 
impurities. 

The degree and nature of the extraction with caustic 
soda determines the alpha cellulose purity in the final 
product, and while there are ways in which alpha cellulose 
may be degraded, it is by caustic soda that it is refined. 
Extraction with caustic soda has itself many variants, and 
plants are operating which carry out this step batchwise, 
continuously, at atmospheric pressure, at superatmospher- 
ic pressure, at low consistency, at high consistency, with 
dilute caustic soda, with strong caustic soda, and with 
combinations of these. The range of temperature em- 
ployed is from 20 to 145°C. 

A simple extraction with hot dilute sodium hydroxide 
will suffice for many purposes. The woodpulp coming 
to this stage comes from a washer-thickener where a large 
part of the water has been removed from it, such that 
it will have three to five parts of water for one of cellu- 
lose. This pulp is then sprayed or mixed with the re- 
quired amount of caustic soda. It may then be mixed 
with steam and enter a continuous tower, where it may be 
held at approximately 95°C. for 2 to + hours, or it may 
be fed to a batch refiner, where it is heated with steam 
at about 20 to 30 pounds pressure (120 to 145°C.) for 15 
to 60 minutes. While these variants are all common, they 
are each important in different ways. It is obvious that 
pressure vessels require less capacity than atmospheric 
pressure vessels, since the time of reaction is shorter, but 
on the other hand they are more expensive to build. Then 
again, some manufacturers add a minimum amount of 
caustic soda and proceed until it is exhausted, whereas 
others operate with an appreciable, constant excess, feeling 
that they gain time and quality uniformity at the expense 
of a little extra chemical. 


In choosing among these processes, an important con- 
sideration is cellulose vield. The work here has shown 
that prolonged caustic extraction can result in serious 
losses in yield with no compensating gain in purity. 
Table 3 and Fi igure 1 show that by prolonging the time 
of extraction at 95°C., a very small adv: antage in caustic 
soda consumption is < achieved, but at a greater sacrifice 
in product yield. 

This phenomenon is very clear at 95°C., and even 
more sharp at lower temperatures, say at 70°C. At higher 
temperatures, W here the reaction is more rapid, the loss 
of yield is not so pronounced when the reactions are kept 
to quite a short time. The proper time for caustic soda 
extraction is a very important matter, on which present 
data are far from sharp, but the present results can be 
summarized as in Table 4+ and Figure 2. 

The Table and Figure suggest that to reach a desired 
level of alpha cellulose, only very short times are required 
at high temperatures, approaching zero time at about 
160°C. This temperature 160°C., is not yet firmly es- 
tablished in the laboratory here. It is known, however, 
that 145° is not enough. It is known also that 170°C. 
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TABLE 3 


YIELD IN Caustic SoDA EXTRACTION 
(95°C., 20% consistency ) 


Alpha ©) NaOH | Time, Yield, Alpha 


Calulose on pulp min, Qi Yield, ‘ 


93 i 60 a 93. 
0 90 Jo 91.5 
5.5 120 5. 91. 


60 os 93.0 
90 ‘ 90.0 
120 Oe 89.5 


60 | o4 90.! 
90 
120 


90 
120 


*Alpt ‘eld Alpha final eld 
- a-yie = yiel¢ 
een Alpha initial — 
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94.5 


60 90 120 
TIME IN MINUTES 
Figure 1—Yield of alpha cellulose in caustic soda extrac- 
tion of chlorinated sulphite wood pulp at 95°C. and 20% 
consistency. (see Table 3). 


carries with it certain dangers, since any residual lignin 
will polymerize at that temperature to form colored 
bodies that are extremely difficult to bleach out. There 
is, however, a temperature at which the time can be ex- 
tremely short, and more extraction simply lowers yield 
with no gain in purity. 

Another variable which should be mentioned is the 
ratio of pulp to water in this stage, the “consistency” 
The evidence suggests that at high consistencies, say 25% 
or over, somewhat higher alpha cellulose is obtained. 
Furthermore, at very high consistencies any given per- 
centage of caustic added to the pulp is at a proportionately 
higher concentration around the fibres. At strengths of 


274. 


TABLE 4 
Caustic SopA EXTRACTION TO 93% ALPHA CELLULOSE 


Temp., Time, NaOH, 
aS min. % on pulp 
60 . ‘. . 
50 6.0 
30 aes 


20 SB. 


6. 





°C 








TEMPERATURE 























TIME IN MINUTES 


Figure 2—Caustie soda extraction of chlorinated sulphite 
wood pulp to 93% alpha cellulose (see Table 4). 


4% caustic soda concentration and over, the fibres them- 
selves are softened, and this softening of the fibres can be 
very important in subsequent use of the cellulose, and is 
essential in making cellulose acetate by most of the com- 
mercial processes in use today. 

Summing up the best conditions of carrying out hot 
caustic extraction: very high consistencies, very short 
times, temperatures about 155 to 165°C., and thorough 
mixing with steam and chemicals. Since the time is very 
short and the consistency high, such equipment should 
be small for a given capacity, and therefore require 
smaller capital outlay. It should be continuous and re- 
quire little if any attention. 

Will this type of refining be successful on the less 
thoroughly digested woodpulp suggested earlier? The 
data show that the degree of purification which can be 
achieved by this procedure is actually greater than can be 
achieved with more thoroughly digested woodpulp. The 
yields based on wood are better. This applies both to 
softwoods and to Northern black cottonwood, a member 
of the genus Populus. 

The degree to which this type of refining (dilute 
caustic soda extraction at high temperatures) is carried 
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out depends on the desired qualities in the finished pro- 
duct. Unfortunately, there is a ceiling of alpha cellulose 
beyond which this type of extraction will not purify cel- 
lulose, doubtless due to competing reactions of purifica- 
tion and degradation'#), To make pure alpha cellulose, 
an extraction with strong caustic soda at low temperature 
is required. The conditions for this extraction have been 
described“), Figure 3 is taken from that paper. It 
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§ 58 2 % 45 55 


TEMPERATURE, °C 


Figure 3—Alpha cellulose obtained by caustic soda ex- 
traction of refined sulphite wood pulp. From Wyman, M., 
and Sherk, D.L., TAPPI, 39, 786 (1956). The lines 
represent conditions for obtaining the levels of alpha 
cellulose represented by the numbers opposite each line. 


can be seen that to achieve 100% alpha cellulose a relative- 
ly narrow set of conditions of temperature and concen- 


tration are required — about 130 g.p.l. NaOH at 1876, 


Just as in hot caustic soda extraction, so also in the 
cold the duration of extraction can seriously affect yield. 
A prolonged extraction can have quite harmful effects. 
Alpha cellulose rises very sharply to a maximum in about 
5 to 10 minutes, then falls off, then climbs again. This is 
an oxidation effect, in which the second inflection re- 
presents the exhaustion of oxygen. It has been found 
possible to eliminate the dip by the use of an organic 
antioxidant. The far cheaper sodium sulphite has not 
been effective here. In addition to alpha cellulose con- 
siderations, the degree of polymerization of the cellulose 
falls off with time in cold caustic extraction, and the 
extraction should therefore be of limited duration. 

This type of cold caustic soda extraction of woodpulp 
has been carried out for a number of years by several 
cellulose manufacturers. The degree of success is gener- 
ally a compromise with caustic soda economy. Wood- 
pulp treated under the desired conditions is swollen and 
slimy and very hard to handle in conventional equip- 
ment. It is very difficult to recover the caustic soda in 
re-useable form. Most systems are at least in part count- 
tercurrent flow(!® !7), involving the formation of a 
pulp sheet from a slurry of highly swollen cellulose. There 
are serious mechanical difficulties associated with achieving 
this at optimal conditions. However, there is equipment 
in operation which treats cellulose with even stronger, 
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caustic soda solution. This is the slurry steeping equip- 
ment widely used in rayon plants throughout the world. 
There are two main embodiments, the Schmitz press (1%) 
being the most common on this continent, and the screw 
press. If pulp wet with about 10% caustic soda is raised 
in consistency to above 40%, (that is 1.5 parts solution 
to 1 part cellulose) the caustic economy becomes tolerable 
without any reuse of caustic soda. 


The type of equipment used in the rayon plants i is too 
small for cellulose plants. What is needed is equipment 
which will treat a cellulose suspension with 10% caustic 
soda solution for 5 to 10 minutes at 25°C., and then thicken 
it to above 40% consistency. This equipment should be 
large enough to fit a modern cellulose plant. It should 
of course be continuous and require little attention. 


Since yields are very much better in cold strong caustic 
soda extraction than they are in hot weak caustic soda 
extraction, the economy of the process is better if the 
main burden of purification is put on the cold system. 
However, there are certain dangers in this. Below a 
certain level of the hot extraction, mannan removal is 
insufficient, and the suitability of the pulp cellulose deri- 
vatives manufacture becomes endangered. This is especi- 
ally so if the original pulping of the wood is carried on 
only to the minimum extent, for highest yield. In 
fact this problem, of removal of mannan, is one of the 
most serious obstacles to the development of a high-yield 
process for acetate-grade cellulose from softwoods (" ), 


Bleaching 


Wood cellulose, after extraction with caustic soda, is 
not white, even after thorough washing, and for nearly 
all uses must be bleached. For some end-product uses, 
especially viscose rayon and cellophane manufacture, con- 
trol of degree of polymerization of the cellulose is also 
highly desirable. These ends are achieved by bleaching 
Ww ith ‘sodium hypochlorite. 


More recently chlorine dioxide has been used success- 
fully in bleaching cellulose, especially where a higher 
degree of polymerization is desired in the final product. 
Again there are different ways of carrying out bleaching, 
some operators preferring a longer time with a minimum 
of bleach—bleaching to exhaustion of bleach—while others 
prefer an excess in a limited time. The best economy and 
control will depend on the equipment available in a given 
plant. 


Deresination 

Wood components which can be very troublesome in 
the finished cellulose are the solvent- soluble extractives 
sometimes collectively called “resin”. Their removal 
from wood cellulose has been discussed recently (2°2!), 
and there is little to add to the previous discussions. The 
use of certain surface-active agents, especially nonyl 
phenol condensed with ethylene oxide, is extraordinarily 
effective when used in the hot caustic soda extraction 
stage for removal of resin. Mixed surface active agents 
are also effective (2). 

Some particular woods give rise to special problems in 
deresination, and various devices have been used to deal 
with them. Earlier it was mentioned that the genus 

Populus had resin particularly difficult to deal with. The 
use of surface active agents will remove the resin, but 
with these woods the quantities required are so large as 
to be prohibitive except for special purposes. However, 


the “resin” in these woods turns out to be highly un- 
saturated and readily susceptible to oxidation. (This 
ready oxidation makes chip storage a serious fire hazard 
Seasoning of the wood has a marked 


with this genus). 
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effect on resin reduction. Another approach to this 
problem is oxidation in place of chlorination as the first 
stage of treatment of the raw pulp, and chlorine dioxide 
has been used for this purpose. It is remarkably effective, 
as shown in Table 5. 


TABLE 5 


EFrect oF ClO. First StaGE TREATMENT ON REMOVAL OF 
RESIN FROM Corronwoop SULPHITE PuLp By Hor 
Caustic Sopa EXTRACTION 
6% NaOH and 0.3% Antarox A-400 on 
cellulose at 140°C. for 20 minutes 


Extractives 
After refining 


First Stage Treatment 
© Chlorine Demand 


70% chlorine......... ‘ 0.467 


35% chlorine + 35% ClO. ‘ 0.364 
60% Cl0..... 0.108 
70% Cl0...... 0.127 
80% Cl0..... 0.118 


Physical form 


In the manufacture of cellulose from wood some 
thought must be devoted to the physical properties of the 
finished product. In most cases cellulose after purifica- 
tion and bleaching is soured with sulphur dioxide solu- 
tion to kill residual bleach and to lower the pH of the 
cellulose in order to release metals. The pH is then 
raiscd, by washing in pure water, or by caustic soda to 
just below neutrality. The pulp is sheeted, and dried, 
and the sheet is cut and packaged for shipment. 

The physical properties required are determined by 
the first stage of treatment in the customers’ plants. In 
cellulose acetate plants woodpulp is usually shredded, and 
it is quite important that the sheet be very ‘soft and readily 
torn or cut apart. For acetate it may not be necessary 
to dry the woodpulp as sheets. There are now several 
machines on the market which will dry the cellulose in 
fluff form, and bale it. This is very similar to cotton, 
and would be quite familiar to the consumers. 


For viscose and most ethers the woodpulp is usually 
steeped in strong caustic soda solution and is required 
to have some rigidity under these conditions. 


If the woodpulp is to be used in slurry steeping, then 
it is important not to have too high a fines content in 
the pulp. For most types of press, the fines tend to go 
with the expressed caustic and cause trouble in the 
recycle or caustic recovery system. There are two types 
of fines, those which are part of the tree as it grows, and 
those which are merely broken fibre fragments. The first 
is very high in resin and ash. The fibre fragments tend 
to be low in alpha cellulose and in degree of polymer- 
ization. The first should be screened out, and the form- 
ation of the second should be avoided by proper choice 
of processing equipment and conditions. Fines repre- 
sent economic loss and quality deterioration. 


In the manufacture of chemical cellulose from wood 
there are today two clear objectives. The first is to make 
a considerable improvement in the economy of cellulose 
manufacture, while maintaining the highest present-day 
standards of purity. Figure 4 shows the present state of 
the art. This paper suggests that this may be improved 
by: 

1. More moderate delignification and defibering of 
the wood, at higher vields, 
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Figure 4—Yield of cellulose during purification 
from coniferous wood. 
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The development of improved purification pro- 
cesses to make the best use of rawer pulp, 


3. Close attention to each stage in processing from 
the point of view of yielé, This may involve the 
dev elopment of new equipment for hot and cold 
caustic soda extraction, 


4. Wider use of members of the genus Populus as a 
source of cellulose. 


The second objective is to make cellulose from wood 
as pure as the best that can be made from cotton. Wood 
cellulose has recently invaded quite successfully areas pre- 
viousy held by cotton ~ cellulose acetate film, triacetate 
yarn and photographic film, molded and extruded plastics 
and others. This has come about because of steady im- 
provement in the quality of wood cellulose. However, 
the best commercial celluloses from wood still are only 
97.5% alpha cellulose and have 0.7 to 0.9% nonglucan 
carbohydrate constituents and 0.05% solvent-soluble ex- 
tractives. By close attention to the conditions necessary 
for removing the last traces of unwanted carbohydrates 
and resins, further advances toward pure cellulose from 
wood will be made. 
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The Settling Behavior of Uranium 


Trioxide- Water Slurries’ 


A. W. BOYD? and J. L. WHITTON? 


Measurements have been made of the settling 
rates of uranium trioxide-water slurries in the range 
of 100 to 800 grams per litre at temperatures of 25 
and 53°C. The applicability of several equations pro- 
posed by others to relate the rate of hindered settling 
to surface area, concentration and viscosity has been 
examined. It is shown that the data are best fitted by 
an empirical equation similar to one proposed by 
Allison and Murray for zinc, aluminum and mag- 
nesium oxide suspensions. 


oe trioxide-water slurries have been proposed 
as fuels for homogeneous nuclear power reactors. The 
settling behavior of such slurries, that is, the variation of 
settling rate with concentration, viscosity and surface area, 
has been investigated because this behavior is of impor- 
tance in the design and operation of slurry reactors. It 
was also thought that settling rate measurements might be 
a convenient way of determining the surface area of the 
uranium trioxide. 

Several equations relating the settling rate of a con- 
centrated suspension to concentration, viscosity and par- 
ticle size have been published. In this paper three of these 
are applied to the results obtained for the uranium tri- 
oxide-water slurries. 


Experimental procedure and results 


Anhydrous uranium trioxide is not stable in water but 
forms one of four hydrates depending on the temperature 
(1), The two that are stable in the temperature range of 
interest in reactors i.e. about 150°-250°C., are the mono- 
hydrates designated as a and , the former consisting of 
rod- -shaped crystals, the latter of square thin crystals. The 
slurries used in the settling rate measurements were pre- 
pared from the two monohydrates of uranium trioxide. 


Pure samples of the a and 6 forms were prepared and 
the particle size distribution of each was measured by a 
combination of sedimentation analysis and microscopic 
analysis. An Andreasen pipette was used in the sedimen- 
tation analysis, the concentration of the suspension being 
1 gm. of oxide per litre and 0.2 gm. of Calgon* (as a 
dispersing agent) per litre. The particle size distribution 


°A trade name for sodium hexametaphosphate. 


1Manuscript received August 20, 1958. 

2Development Chemistry Branch, Atomic Energy of Canada Limited, 
Chalk River, Ont. 

Based on a paper presented at the 40th Annual Conference, The Chemical 
Institute of Canada, Vancouver, B.C., June 1957. 
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curve thus obtained for the B-UO,.H,O is shown in 
Figure 1. In order to obtain the distribution in the region 
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Figure 1—Sedimentation analysis of G-UO;.H.O in water. 


below 10% of the weight undersize, the particles in the 
range below 4 microns were measured using a microscope 
with a filar micrometer. About 500 particles were 
measured in each case. In order to obtain the surface area 
the two measurements were combined and it was calcu- 
lated that the surface area of the a-UO,.H,O was 1.5 x 
10* cm.? per cm.’ and that of the £- UO, H, ,O was 6.1 x 
103 cm.? per cm.*. These values correspond. to rods (a) 
averaging 10 x 1.5 x 1.5 microns and to plates (f) averag- 
ing 13 x 13 x 3 microns. 


The settling measurements were made using graduated 
glass cylinders (i.d. 2.4 cm.) immersed in a thermostated 
water bath. Measurements were made at 25°C. and 53°C. 
using several concentrations of each of the monohydrates. 
The lowest concentration used was the minimum needed 
to obtain a sharp interface between the sediment and the 
supernatant. Mixing was done by rotation of the graduates 
and all measurements were repeated. 


The amount of uranyl nitrate in solution in the slurries 
used was of the order of 10 parts per million. The value 
used for the density of both hydrates was 6 gm./cm.3, 


277 































30 


volume (mt.) 





20 


Settled 


400g/L 


- (mi.) 





Volume 


20 


tried 





300 g/.. 





Figure 2—-Settling curves of G-UO;.H.O in water at 25°C. 
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Figure 4—-Settling curves of a-UO;.H.O in water at 25°C. 


The results of the settling rate measurements are shown 
in Figures 2, 3, 4 and 5, where the settled volume is 
plotted against time. 

The shapes of the settling curves shown in Figures 2 to 
5 are of interest as they illustrate two distinct. types of 
settling. In the first case, where the rate of settling remains 
constant over a wide range of concentration, the increase 
in concentration must be occurring at the bottom of the 
suspension and is propagated upward at a rate that is slow 
compared to the settling rate. Hence the top of the 
suspension remains at the initial concentration. In the 
second case, the rate of settling decreases with increasing 
concentration and here the concentration throughout the 
whole height of the suspension must be changing as it 
settles. This is shown in the curve for the 300 gm./litre 
slurry in Figure 2. 

While the initial concentration determines the velocity 
at any point in the first type of settling, it is not so 
obvious what the effective concentration at any point is 
when the settling curve is of the second type. An equation 
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Figure 3—Settling curves of G-UO;.H:0 
in water at 53°C. 








( seconds) 


Figure 5—Settling curves of a-UO;.H.O in water at 53°C. 


relating effective slurry concentration to the shape of the 
settling curve has been derived by Kynch ©) in which it 
is shown that the ratio of the height of the settled 
suspension to the time of settling is directly related to the 
concentration. Consider the curve shown in Fi igure 6. The 
effective concentration, that is the concentration that 
determines the settling rate, at the point P is the con- 
centration corresponding to the intercept on the Y axis 
of the slope at P. Thus in Figure 6 the effective concentra- 
tion at P is twice the initial concentration of the sus- 
pension. Kynch’s equation can be written in the form 


where T is the height of the intercept on the Y axis and 
x and t are the coordinates of the point P. In order to 
determine how well the measurements of the settling 
rates of the slurries conformed to Kynch’s equation the 
following experiment was carried out. 
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Figure 6—Kynch’s method of determining effective con- 
centration. 


A slurry with a a settling curve similar to that of the 
300 poy A slurry in Figure 2 was used. Tangents to 
this curve were drawn from various intercepts on the Y 
axis and the slopes measured. Slurries with concentrations 
corresponding to the various intercepts were then pre- 
pared and their initial settling rates determined. A com- 
parison of the values of the slopes of the tangents and the 
values of the initial settling rates of the slurries is shown 
in Table 1. 


TABLE 1 


COMPARISON OF OBSERVED AND CALCULATED* SETTLING RATES 


Observed eres Cale ulated Settling 
Rate cm./sec. x 104 Rate cm./sec. x 104 


Concentration 
gm./litre 


617 373 345 
800 | 215 | 214 
1000 118 | 121 
1300 | 56 56 
1600 33 | 32 


*Calculated from Kynch’s equation. 


These data show that the measured settling rate of a 
suspension of this type may be used to predict accurately 
the settling rates of various concentrations of the same 
suspension. 


Application of the Steinour equation 


The following equation has been obtained empirically 
by Steinour ‘) and is equivalent at low concentrations to 
an equation derived by Hawksley ‘*): 


V eV 2 1088-0. en, (2) 


This equation has recently been applied to the settling of 
thorium oxide-water slurries by Reed and Crowley ) 
who plotted their data as log V against log <* 10°1.83(1-€), 
If the volume of the solid in the slurry is increased by the 
absorption of water, the effective volume of the solid is 
(1—e€)(1+ a) where a is the ratio of the volume of 
the liquid absorbed to the volume of solid), but the log-log 
plot should still result in a straight line. 
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Figure 7—Measured settling rate plotted against e? 
19-'-820'-© (Steinours equation ) 


- 


In Figure 7 the settling rates derived from Figures 2 
to 5 are plotted against <? 10-!-82(1-© on a log-log scale. 
When ¢ is equal to unity, V should equal V, and accord- 
ing to Figure 7 the Stokes velocity should be greater for 
the a-UO,.H,O than for the 6- UO, .H,O. The ratio of 
the measured values of the surface areas of the a-UO,.H,O 
and §-UO,.H,O is 2.5:1 however, indicating that the 
settling rate of the a- compound should be about one 
sixth that of the B-compound. 


The Steinour equation fails to fit the data because it 
does not give a rapid enough decrease in settling rate with 
increasing concentration. 


Application of Loeffler and Ruth equation 


The Loeffler and Ruth ‘®) equation for the settling of 
suspensions is based on the Kozeny equation for the flow 
of fluid through a porous bed, w hich is as follows: 


g é 
V= MER ies a ace ite tae 3 
5n S?(1 — €)? “ 


In applying this to suspensions, Loeffler and Ruth substi- 
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Figure 8 -Measured settling rate plotted against 
e*/1 t 
11 + €/1 € 
(Loefiier and Ruth’s equation) | 





tute the difference in density between solid and liquid 
minus the viscous drag on the particles for AP and get: 





where V, is the Stokes velocity and k and K depend on 
the ” ape of the particles in the slurry. 

n Figure 8 the settling rates derived from Figures 2 
to § are plotted against 















Phe value of 11 for k in the denominator was obtained 
from the observed shapes of the particles. If the LoefHler 
and Ruth equation were 
Figure 
origin. 


applicable here, che points in 
should lie on straight lines that go though the 


This equation, like the Steinour one, does not fit the 
data because it also does not give a rapid enough decrease 
in settling rate with increasing concentration. 


Application of the Allison and Murray equation 


Allison and Murray “ start with the Kozeny equation 
and substitute d,, the density of the settling suspension, 
for AP, the pressure drop across the bed. While the 
correct substitution must include the density difference 
between the solid and the liquid, the equation they obtain 


v € 


5n S? (1 € 


fits their data and the data derived from Figures 2 to 5 


quite accurately. In Figure 9 the rates of settling from 


> 


Figures 2 to 5 are plotted against The term d, 


\ ¢}* 
omitted as it does not appear to have any theoretical basis 
and a better fit is obtained without it. 
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Figure 9 — Measured settling rate plotted against €°/(1 — ¢*) 
modified Allison and Murray equation. ) 


Figure 9 shows that the dependence of settling rate 
on concentration is given quite accurately by the function 

(1—e)*. In order to determine if the settling rate is 
inversely proportional to the viscosity and the square of 
the surface area, the ratios of the slopes of the lines in 
Figure 9 are compared with the ratios of measured surface 
areas and viscosities in Table 2. 


ABLE 2 


COMPARISONS OF MEASURED AND CALCULATED 


€ 3 


Ratios of Kin V = K 


(1 


Calculated Measured 






Ks 

: 1.71 1.49 B-UO;.H:0 
K 259 1.68 a-UO;.H,0 
. 

a 6.0 8.8 25°C. 
Ke 10.0 53°C. 


For the uranium trioxide slurries, this equation at least 
gives a reasonable indication of the effect of changes in 
viscosity and surface area on settling rate. Allison and 
Murray’s measurements were made using suspensions of 
aluminum, zinc and magnesium oxides in water. The 
range of porosities was approximately the same as in the 
uranium trioxide-water slurries and the size of the par- 
ticles in these oxides was of the order of a few microns. 





Conclusion 


Uhe data show that for uranium trioxide-water slurries, 
the dependence of settling rate on concentration is 
«'/(1—e)* and suggest that this equation may also apply 
to other suspensions consisting of non-spherical micron- 
sized particles over a porosity range of 0.85-0.95. 
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